Elementary Number Theory 


A revision by Jim Hefferon, St Michael’s College, 2003-Dec 
of notes by W. Edwin Clark, University of South Florida, 2002-Dec 



source compiled on January 5, 2004 by Jim Hefferon, jim@joshua.smcvt.edu. 


License restriction claimed by W. Edwin Clark. Copyleft 2002: “ Copyleft means that 
unrestricted redistribution and modification are permitted, provided that all copies 
and derivatives retain the same permissions. Specifically no commerical use of these 
notes or any revisions thereof is permitted.” 





Preface 


Mathematics is the queen of sciences and 
arithmetic the queen of mathematics 

Carl Friedrich Gauss 

Number theory, known to Gauss as “arithmetic,” studies the properties of the 
integers: ... — 3, —2, —1, 0,1,2,3 .... Although the integers are familiar, and 
their properties might therefore seem simple, it is instead a very deep subject. 

For example, here are some problems in number theory that remain unsolved. 
(Recall that a prime number is an integer greater than 1 whose only positive 
factors are 1 and the number itself.) Note that these problems are simple to 
state just because a topic is accessibile does not mean that it is easy. 

1. ( Goldbach’s Conjecture) Is every even integer n > 2 the sum of two primes? 

2. ( Twin Prime Conjecture) Are there are infinitely many twin primes? 
( Twin primes differ by 2, like 11 and 13.) 

3. Are there infinitely many primes of the form n 2 + 1? Of the form 2" — 1? 
(Ones of this form are Mersenne primes.) Of the form 2 2 + 1? (These 
are Fermat primes.) 

4. Are there infinitely many primes whose digits are all l’s? (Numbers of 
this form are repunits.) 

5. Are there infinitely many perfect numbers? (An integer is perfect if it is 
the sum of its proper divisors; 6 is perfect because 1 + 2 + 3 = 6.) 

6. (3 n + 1 Conjecture) Consider the function / defined by: f(n) = 3n + 1 
if n is odd and f(n) = n/2 if n is even. Does the sequence of iterates 
/(n), /(/(n)), /(/(/(«))),... always contain 1, no matter what starting 
value n is used? 

7. Is there a fast algorithm for factoring large integers? 

So the subject is not simple. However it is accessible, and beautiful. 
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PREFACE 


A caution In some areas a person needs to learn by starting from first princi¬ 
ples. The first course in Calculus is like that; students learn limits first to avoid 
getting nutty ideas about nx n ~ , But other areas are best mastered by diving 
right in. 

In this book you dive into mathematical arguments. Number Theory is right 
for this in part because of its accessibility. 

But always keep in mind the caution: do not underestimate the material. 
You will find this subject hard, albiet rewarding. 

Prerequisites We require only Calculus I. Even that requirement is not strict 
(limits come up, as do the rules of logarithm manipultion), so the main purpose 
of the prerequisite is that we expect that with it comes a certain amount of 
mathematical maturity, including familiarity with basic set theory and some 
function facts. 

Other resources The Internet contains much interesting and current infor¬ 
mation about number theory; see the Bibliography. The websites by Chris 
Caldwell [2] and by Eric Weisstein [13] are especially good. To see what is going 
on at the frontier of the subject, you may take a look at some recent issues of 
the Journal of Number Theory which you will find in any university library. 
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Chapter 1 

Divisibility 


In this book, all numbers are integers, unless specified otherwise. Thus in the 
next definition, d, n, and k are integers. 

1.1 Definition The number d divides the number n if there is a k such that 
n = dk. (Alternate terms are: d is a divisor of n, or d is a factor of n, or n is 
a multiple of d.) This relationship between d and n is symbolized d \ n. The 
symbol d j n means that d does not divide n. 

Note that the symbol d \ n is different from the fraction symbol d/n. It 
is also different from n/d because d \ n is either true or false, while n/d is a 
rational number. 

1.2 Theorem (Divisibility Properties) For all numbers n, m, and d, 

(1) d | 0 

(2) 0 | n => n = 0 

(3) 1 | n 

(4) (Reflexivity property) n \ n 

(5) n | 1 => n = 1 or n = — 1 

(6) (Transitivity property) d \ n and n \ in => d \ m 

(7) (Multiplication property) d \ n =t- ad \ an 

(8) (Cancellation property) ad \ an and a/0 => d \ n 

(9) (Linearity property) d \ n and d \ m => d \ an + bm for all a and b 

(10) (Comparison property) If d and n are positive and d \ n then d < n 

PROOF. For the first item, take k = 0. For the second, if 0 | n then n = 0 • k = 0. 

The next item holds because we can take n as the k in the definition. Re¬ 
flexivity is similar: n = n ■ 1 shows that it holds. The next property follows 
immediately from Basic Axiom 3 for Z, from the first Appendix. 

For Transitivity, assume the d \ n and that n \ m. Then n = dk\ and 
to = nfc 2 for some ki,k% £ Z. Substitute to get to = nfe = (dk\)k 2 - By the 
Associative Property of Multiplication, (dk i)/c 2 = ^(fcife), which shows that d 
divides to. 
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Multiplication also follows from associativity. Assume that d \ n so that 
n = dk. Then an = a(dk) = ( ad)k shows that ad \ ak. 

For Cancellation, assume that a yf 0 and that ad \ an. Then there is a k 
such that an = ( ad)k. We will show that n = dk. Assume first that a > 0. By 
the Trichotomy Property from the first Appendix, either n > dk or n = dk or 
n < dk. If n > dk then we have that an > a(dk) = ( ad)k , which contradicts this 
paragraph’s assumption that an = ( ad)k. If n < dk then an < a(dk) = ( ad)k, 
also contradicting the assumption. Therefore n = dk, and so d j n. The 
argument for the a < 0 case is similar. 

To verify Linearity, suppose that d \ n and d \ m so that n = dk\ and 
to = dk2 for k\, /c 2 € Z. Then an + bm = a(dki) + b(dk2) = d(ak\ + bk2) shows 
that d | (an + bm). 

Finally, for Comparison, assume that d, n > 0 and d \ n. Then n = dk 
for some k. Observe that k is positive because the other two are positive. By 
Trichotomy, either d < n or d = n or d > n. We will show that the d > n case 
is not possible. Assume that d > n. Then dk > nk follows by one of the first 
Appendix’s Properties of Inequalities. But that gives n > nk, which means that 
n ■ 1 > n ■ k despite that fact that k is positive and so 1 < k. This is impossible 
because it violates the same Property of Inequalities. QED 

1.3 Definition An integer n is even (or has even parity) if it is divisible by 2 
and is odd (or is of odd parity) otherwise. 

1.4 Lemma Recall that |a| equals a if a > 0 and equals —a if a < 0. 

(1) If d | a then —d \ a and d \ —a. 

(2) If d | a then d \ |a| 

(3) The largest positive integer that divides a nonzero number a is |a|. 

PROOF. For (1), if d | a then a = dk for some k. It follows that a = (— d)(—k ) 
and since — d and —k are also integers, this shows that — d \ a. It also follows 
that —a = (— k)d, and so d \ —a. 

For (2), suppose first that a is nonnegative. Then |a| = a and so if d \ a 
then d \ \a\. Next suppose that a is negative. Since |a| = —a for negative a, and 
since (1) shows that d \ —a, and d therefore divides |a|. 

For (3), first note that |a| actually divides a: in the a > 0 case |a| | a because 
in this case |a| — a and we know that a \ a, while in the a < 0 case we have that 
a = |a|(—1), so that |a| is indeed a factor of a. We finish by showing that |a| 
is maximal among the divisors of a. Suppose that d is a positive number that 
divides a. Then a = dk for some k, and also —a = d(—k). Thus d \ |o|, whether 
a is positive or negative. So by the Comparison property of Theorem 1.2, we 
have that d < |a|. qed 



Chapter 2 

Prime Numbers 


2.1 Definition An integer p > 2 is prime if it has no positive divisors other 
than 1 and itself. An integer greater than or equal to 2 that is not prime is 
composite. 

Note that 1 is neither prime nor composite. 

2.2 Lemma An integer n > 2 is composite if and only if it has factors a and 
b such that 1 < a < n and 1 < b < n. 

PROOF. Let n > 2. The ‘if’ direction is obvious. For ‘only if’, assume that n is 
composite. Then it has a positive integer factor a such that a ^ 1, a ^ n. This 
means that there is a b with n = ab. Since n and a are positive, so is b. Hence 
1 < a and 1 < b. By Theorem 1.2, a < n and b < n. Since a/1 and a ^ n we 
have 1 < a < n. If b = 1 then a = n, which is not possible, so b ^ 1. If b = n 
then o=l, which is also not possible. So 1 < b < n, finishing this half of the 
argument. QED 

2.3 Lemma If n > 1 then there is a prime p such that p \ n. 

PROOF. Let S denote the set of all integers greater than 1 that have no prime 
divisor. We must show that S is empty. 

If S is not empty then by the Well-Ordering Property it has a smallest 
member; call it m. Now m > 1 and has no prime divisor. Then m cannot be 
prime (as every number is a divisor of itself). Hence m is composite. Therefore 
by Lemma 2.2, m = ab where 1 < a < m and 1 < b < m. Since 1 < a < m, the 
factor a is not a member of S. So a must have a prime divisor p. Then p \ a 
and a | m, so by Theorem 1.2, p \ to. This contradicts the assumption that ?n 
has no prime divisor. So the set S must be empty. QED 

2.4 Theorem (Euclid’s Theorem) There are infinitely many primes. 

PROOF. Assume, to get a contradiction, that there are only a finitely many 
primes pi = 2, p 2 = 3,..., p n . Consider the number N = P 1 P 2 • ■ - p n + 1. 
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Since p± > 2, clearly N > 2. So by Lemma 2.3, N has a prime divisor 
p. That prime must be one of pi , ... , p n since that list was assumed to be 
exhaustive. However, observe that the equation 

N = pi (jpip 2 • • -Pi-iPi+i ■■■p n ) + 1 

along with 0 < 1 < Pi shows by Lemma 3.2 that n is not divisible by p.j. This is 
a contradiction; it follows that the assumption that there are only finitely many 
primes is not true. QED 

2.5 Remark Eucild’s Theorem, and its proof, is often cited as an example of 
the beauty of Mathematics. 

2.6 Theorem If n > 1 is composite then n has a prime divisor p < yfn. 

PROOF. Let n > 1 be composite. Then n = ab where 1 < a < n and 1 < b < n. 
We claim that at least one of a or b is less than or equal to yfn. For if not then 
a > yfn and b > y/n, and hence n = ab > yfn ■ yfn = n, which is impossible. 

Suppose, without loss of generality, that a < yfn.. Since 1 < a, by Lemma 2.3 
there is a prime p such that p \ a. Hence, by Transitivity in Theorem 1.2, since 
a | n we have p | n. By Comparison in Theorem 1.2, since p | a we have 
p < a < yfn. QED 

We can use Theorem 2.6 to help compute whether an integer is prime. Given 
n > 1, we need only try to divide it by all primes p < yfn. If none of these 
divides n then n must be prime. 

2.7 Example Consider the number 97. Note that \/97 < -\/K)0 = 10. The 
primes less than 10 are 2, 3, 5, and 7. None of these divides 97, and so 97 is 
prime. 



Chapter 3 

Division 


3.1 Theorem Where a and b > 0 are integers, there are integers q and r, 
called the quotient and the remainder on division of a by 6, satisfying these two 
conditions. 

a = bq + r 0 < r < b 
Further, those integers are unique. 

Note that this result has two parts. One part is that the theorem says there 
exists a quotient and remainder satisfying the conditions. The second part is 
that the quotient, remainder pair are unique: no other pair of numbers satisfies 
those conditions. 


PROOF. To verify that for any a and b > 0 there exists an appropriate quotient 
and remainder we need only produce suitable numbers. Consider these. 


q = 



r = a — bq 


Obviously a = bq + r, so these satisfy the first condition. To finish the existence 
half of this proof, we need only check that 0 < r < b. The Floor Lemma from 
the Some Properties of R appendix gives 


a 

b 


- 1 < 



< 


a 

b' 


Multiply all of the terms of this inequality by —b. Since b is positive, —b is 
negative, and so the direction of the inequality is reversed. 


b — a > —b 



> —a 


Add a to all three terms of the inequality and replace [a/b\ by q to get 


b > a — bq > 0. 


Since r = a — bq this shows that 0 < r < b. 
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We still must prove that q and r are unique. Assume that there are two 
quotient, remainder pairs 


a = bqi + r x 

with 

0 < 7*1 < b 

a = bq 2 + r 2 

with 

0 < r 2 < b. 


Subtracting 

0 = a - a = {bq 1 + n) - ( bq 2 + r 2 ) = 6 (gi - q 2 ) + 0 i - r 2 ) 
implies that 

(3.1) r 2 - ri = 6 (gi - q 2 ). 

We must show that the two pairs are equal, that r\ = r 2 and qi = q 2 . To obtain 
a contradiction, suppose otherwise. First suppose that ri ^ r 2 . Then one must 
be larger than the other; without loss of generality assume that r 2 > r\. Then 

0 < j*i < r 2 < b 

and so r 2 — r\ < b. But (3.1) shows that b divides r 2 — r\ and by the Comparison 
property of Theorem 1.2 this implies that b < r 2 — r\. This is the desired 
contradiction and so we conclude that 7*1 = r 2 . With that, from equation 3.1 we 
have 0 = b{q\ — q 2 ). Since b > 0, this gives that q\ — q 2 =0 and so q\ = q 2 . QED 

3.2 Corollary The number d divides the number n if and only if on division 
of n by d the remainder is 0 . 

PROOF. If the remainder is 0 then n = dq + 0 = dq shows that d \ n. For the 
other half, if d \ n then for some k we have n = dk = dk + 0 (with 0 < 0 < d) 
and the fact that the quotient, remainder pair is uniqus shows that k and 0 
must be the quotient and the remainder. QED 

That corollary says that Theorem 3.1 generalizes the results on divisibility. 
For instance, fix b = 3. Then, given a, instead of only being able to say that 
a is divisible or not, we can give a finer description: a leaves a remainder of 0 
(this is the case where b \ a), or 1 , or 2 . 

3.3 Definition For b > 0 define a mod b = r where r is the remainder when a 
is divided by b. 

For example 23 mod 7 = 2 since 23 = 7-3 + 2 and —4 mod 5 = 1 since 
-4 = 5- (—1) + 1. 




Chapter 4 


Greatest Common Divisor 


4.1 Definition An integer is a common divisor of two others if it divides both 
of them. 

We write C(a, 6) for the set of numbers that are common divisors of a and b. 

4.2 Definition The greatest common divisor of two nonzero integers a and 6, 
gcd(a, b), is the largest integer that divides both, except that gcd(0, 0) = 0. 

The exception is there because every number divides zero, and so we specially 
define gcd(0,0) to be a convienent value. 

4.3 Example The set of common divisors of 18 and 30 is 

C(18,30) = {-1,1,-2,2,-3,3,-6,6}. 


So, gcd(18,30) = 6. 

4.4 Lemma gcd(a, 6) = gcd(6, a). 

PROOF. Clearly the two sets C{a,b ) and C{b,a) are equal. It follows that their 
largest elements are equal, that is, that gcd(a, b) = gcd(6, a). QED 

4.5 Lemma gcd(a, b) = gcd(|a|, |6|). 

PROOF. If a = 0 and 6=0 then |a| = a and |6| = 6, and so in this case 
gcd(a, 6) = gcd(|a|, |6|). Suppose that one of a or 6 is not 0. Lemma 1.4 shows 
that d | a <t=> d | |a|. It follows that the two sets C(a,b ) and C(\a\, |6|) are the 
same set. So the largest member of that set, the greatest common divisor of a 
and 6 , is also the greatest common divisor of |a| and | 6 |. QED 

4.6 Lemma If a ^ 0 or 6 ^ 0, then gcd(a, 6) exists and satisfies 

0 < gcd(a, 6) < min{|o|, |6|}. 
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PROOF. Note that gcd(a, b) is the largest integer in the set C(a,b). Since 1 | a 
and 1 | b we know that 1 £C(a,b). So the greatest common divisor must be at 
least 1, and is therefore positive. On the other hand, if d € C(a,b) then d | |a| 
and d \ \b\, so d is no larger than |a| and no larger than \b\. Thus, d is at most 
the minimum of \a\ and |6|. QED 

4.7 Example The above results give that 

gcd(48,732) = gcd(—48, 732) = gcd(-48, -732) = gcd(48, -732). 

We also know that 0 < gcd(48,732) < 48. Since if d = gcd(48, 732) then d | 48, 
to find d we need check only for positive divisors of 48 that also divide 732. 

4.8 Remark Observe that the first two lemmas, which draw conclusions about 
the properties of the gcd operator, preceed Lemma 4.6, which shows that the 
gcd exists. 

If two numbers have a greatest common divisor of 1 then they have no 
nontivial common factors. 

4.9 Definition Two numbers are relatively prime if they have a greatest com¬ 
mon divisor of 1. 

Although the relatively prime relationship is symmetric — if gcd(a, b) = 1 
then gcd(6, a) = 1 — we sometimes state it as “a is relatively prime to 6.” 

4.10 Lemma If g = gcd(a, b) then gcd(a/g, b/g) = 1. 

PROOF. The greatest common divisor of a/g and b/g must exist, by the prior 
result. Let gcd(a/g, b/g) = k. Then k is a divisor of both a/g and b/g so there 
are numbers j a and jb such that j a k = a/g and jbk = b/g. Therefore j a {kg ) = a 
and jb(kg) = b , and so kg is a common divisor of a and b. If k > 1 this would 
be a contradiction, because then kg > g but g is the greatest common divisor. 
Therefore k = 1. QED 



Chapter 5 

Bezout’s Lemma 


5.1 Definition A number c is a linear combination of the numbers a and b if 
c = as + bt for some s and t. 

5.2 Example The number c = 16 is a linear combination of a = 2 and b = 5 
since taking s = 3 and t = 2 gives 16 = 3 • 2 + 2 • 5. The number 21 is not a 
linear combination of 2 and 4, since in the equation 21 = s ■ 2 + t ■ 4 the right 
side is divisible by 2 while the left is not. (That is, there are no integers s and 
t; we can solve the equation with rational numbers that are not integral.) 

Thus, the Linearity statement in the Divisibility Properties theorem says 
that if d divides a and 6, then d divides all linear combinations of a and b. So 
the greatest common divisor, gcd(a, b), divides every member of £(a, b ), the set 
of all linear combinations of a and b. The next result says that gcd(a, b) is itself 
a member of C(a, b). 

5.3 Lemma (Bezout’s Lemma) The greatest common divisor of two num¬ 
bers is a linear combination of those two: for all integers a and b there exist 
integers s and t such that gcd(a, b) = sa + tb. 

PROOF. If a and b are 0 then s and t may be anything since gcd(0,0) = 0 = 
s ■ 0 + t ■ 0. So we may assume that « / 0 or fi ^ 0. Consider the set C(a, b) = 
{na + mb : n,m £ Z} of all linear combinations of a and b. 

Denote the set of positive members of £(a, b) by £ + (a, b). Note that £(a, b) 
contains a , —a, b and —b, and since a ^ 0 or b ^ 0, at least one of these four 
numbers is positive. Therefore C + (a, b) is not empty. Because of this, by the 
Well-Ordering Property for N, we know that the set £ + (a, b) contains a smallest 
positive integer; call it d. We will show that d is the greatest common divisor 
of a and b. That will finish the argument because d is a linear combination of 
a and b as it is a member of C. 

Since d £ £ + (a,b) we have d = sa + tb for some integers s and t. Let 
g = gcd(a, b). Then g \ a and g | 6, so by the Linearity property of Theorem 1.2, 
we have g | ( sa+tb ), that is, g \ d. Since g and d are positive, by the Comparision 
property of Theorem 1.2, we have that g < d. 
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If we show that d is a common divisor of a and b, then we will have that 
d < g (as g is the greatest of the common divisors), and so we will have shown 
that g = d. To show that d \ a, write a = dq + r where 0 < r < d and compute 

r = a — dq — a — (sa + tb)q = (1 — sq)a + (— tq)b. 

to conclude that r £ £(a, b). Thus, if r were to be strictly greater than 0 then r 
would be a member of C + (a, b). But this cannot be, since r is strictly less than 
d and d is the smallest integer in C + (a, b). So we must have that r = 0. That 
is, a — dq , and hence d \ a. A similar argument shows that d \ b. Thus, d is 
indeed a common divisor of a and b , and d = g = gcd(a, b). QED 

5.4 Example Notice that 1 = gcd(2, 3) and 1 = (—1)2 + 1-3. Notice also 
that 1 = 2-2+ (—1)3. So the numbers s and t in Bezout’s Lemma are uniquely 
determined. In fact, as we will see later that for each pair a, b there are infinitely 
many s and t. 

5.5 Corollary The set C(a. b) of all linear combinations of a and b equals the 
set of multiples of gcd(a, b). 

PROOF. We observed above that any member of C(a, b ) is a multiple of gcd(a, b ). 
For the converse, consider the multiple fc-gcd(a, 6), apply Bezout’s Lemma to get 
s, t € Z so that gcd(a, b) = sa + tb, and substitute: k ■ gcd(a, b) = k ■ (sa + tb) = 
(ks)a+(kt)b. QED 

5.6 Lemma If a | be and a is relatively prime to b then a \ c. 

PROOF. Since gcd(a, b) = 1, by Bezout’s Lemma there are coefficients s and t 
such that 1 = as+bt. Multiply both sides by c to get c = cas+cbt = a(cs)+(bc)t. 
Note that a \ a(cs) and that a | be by assumption, so Theorem 1.2 gives that a 
divides the linear combination a(cs) + ( bc)t = c. QED 

Observe that 6 | (4-9) but 6 j 4 and 6 f 9 (6 is not relatively prime to 4, 
and is also not relatively prime to 9). Thus the condition of relative primality 
is needed in that lemma. 

We can completely characterize C(a,b). 

5.7 Lemma Fix a,b £ Z. If gcd(a, b) \ c then the equation sa + tb = c has 
infinitely many solution pairs s, t , which have the form 

s=s 0 -j ■ ( b/d ), t=t 0 +j- ( a/d) j £ Z 

where sqAo is any particular solution pair. 

PROOF. First assume that a solution pair soAo exists, to show that any pair of 
numbers of that form also solve the equation. Plug them into the equation. 

(so-j(b/d))-a+(t 0 +j(a/d))-b = ( s 0 a+t 0 b)+j(-(ab/d)+(ab/d )) = s 0 a+t 0 b = c 
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To finish we must show that pairs of the stated type are the only solutions.. 
Suppose that s and t also solve the equation: sa + tb = c. Subtracting gives 
(s — so)a + (t — to)b = 0, that is, 

(*) (s - s 0 )a = (t 0 - t)b. 

Divide by g = gcd(a, b) on both sides to get (s — s 0 )(a/g) = (to — t)(b/g ), which 
shows that b/g divides (s — so)(a/g). By Lemma 4.10, gcd (a/g,b/g) = 1 and 
so the prior result, Lemma 5.6, shows that b/g divides s — So- Thus, for this 
solution pair s, t, there is a j £ Z such that j ■ (b/g) = s — sq, that is, s has the 
form s = So ~ j (b/g)- With that form for s, substituting into equation (*) gives 
that ((so — j(b/g)) — s 0 )a = — ja(b/g ) equals (t 0 — t)b. Dividing both sides by 
b and rearranging gives that t = t 0 + j (a/g). QED 
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Chapter 6 

The Euclidean Algorithm 


We can efficiently compute the greatest common divisor of two numbers. 

We first reduce the problem. Since gcd(a,6) = gcd(|a|, |fo|) (and gcd(0,0) = 
0), we need only give a method to compute gcd(a, b) where a and b are nonneg¬ 
ative. And, since gcd(a, b) = gcd(6, a), it is enough for us to give a method for 
a > b > 0. 

6.1 Lemma If a > 0 then gcd(a, 0) = a. 

PROOF. Since every integer divides 0, C(a, 0) is just the set of divisors of a. The 
largest divisor of a is |a|. Since a is positive, |a| = a, and so gcd(a,0) = a. QED 

The prior lemma reduces the problem of computing gcd(a, b) to the case 
where a > b > 0. 

6.2 Lemma If a > 0 then gcd(a, a) = a. 

PROOF. Obviously, a is a common divisor. By Lemma 4.6, gcd(a, a) < \a\ and 
since a is positive, |a| = a. So a is the greatest common divisor. QED 

We have now reduced the problem to the case a > b > 0. The central result 
is next. 

6.3 Lemma Let a > b > 0. If a = bq + r, then gcd(a, b) = gcd(6, r). 

proof. It suffices to show that the two sets C(a, b) and C(b, r) are equal, because 
then they must have the same greatest member. To show that the sets are equal 
we will show that they have the same members. 

First, suppose that d € C(a,b), so that d \ a and d \ b. Note that r = a — bq. 
By Theorem 1.2(3) we have that d \ r. Thus d \ b and d \ r, and so d £ C(b,r). 
We have shown that any member of C(a, b) is a member of C(5, r), that is, that 
C(a , b) C C(&, r). 

For the other containment, assume that d € C(b,r) so that d \ b and d \ r. 
Since a = bq + r, Theorem 1.2(3) applies again to shows that d\ a. So d \ a and 
d | 6, and therefore d € C(a,b). QED 
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CHAPTER 6. THE EUCLIDEAN ALGORITHM 


The Euclidean Algorithm uses Lemma 6.3 to compute the greatest common 
divisor of two numbers. Rather introduce a computer language in which to give 
algorithm, we will illustrate it with an example. 


6.4 Example Compute gcd(803,154). 


gcd(803,154) = gcd(154,33) 
gcd(154,33) = gcd(33, 22) 
gcd(33,22) = gcd(22,11) 
gcd(22,11) = gcd(ll, 0) 
gcd(ll, 0) = 11 


since 803 = 154-5 + 33 
since 154 = 33 • 4 + 22 
since 33 = 22 • 1 + 11 
since 22 = 11 • 1 + 0 


Hence gcd(803,154) = 11. 

6.5 Remark This method is much faster than finding C(a,b) and can find 
gcd’s of quite large numbers. 

Recall that Bezout’s Lemma asserts that given a and b there exists two 
numbers s and t such that gcd(a, b ) = s-a + t-b. We can use Euclid’s Algorithm 
to find s and t by tracing through the steps, in reverse. 

6.6 Example Express gcd(803,154) as a linear combination of 803 and 154. 

11 = 33 + 22- (-1) 

= 33 + (154 - 33 • 4) • (-1) = 154 • (-1) + 33-5 
= 154 • (-1) + (803 - 154 • 5) • 5 = 803 • 5 + 154 • (-26) 



Chapter 7 

The Fundamental Theorem 


7.1 Theorem (Fundamental Theorem of Arithmetic) Every number 
greater than 1 factors into a product of primes n = P 1 P 2 • • -p s - Further, writ¬ 
ing the primes in ascending order pi < p 2 < • • • < p s makes the factorization 
unique. 

Some of the primes in the product may be equal. For instance, 12 = 2-2-3 = 
2 2 -3. So the Fundamental Theorem is sometimes stated as: every number greater 
than 1 can be factored uniquely as a product of powers of primes. 

7.2 Example 600 = 2 • 2 • 2 • 3 • 5 • 5 = 2 3 • 3 • 5 2 

We will break the proof of the Fundamental Theorem into a sequence of 
Lemmas. 

7.3 Lemma (Euclid’s Lemma) If p is a prime and p \ ab, then p \ a or p \ b. 

PROOF. Assume that p \ ab. If p \ a then we are done, so suppose that it does 
not. Let d = gcd(p, a). Note that d > 0, and that d \ p and d \ a. Since d \ p 
we have that d = 1 or d = p. If d = p then p \ a, which we assumed was not 
true. So we must have d = 1. Hence gcd(p, a) = 1 and p \ ab. So by Lemma 5.6, 
p | b. QED 

7.4 Lemma Let p be prime. Let a\, a 2 ,..., a n , n > 1, be integers. If p \ 
aia 2 • • • a n , then p \ ai for at least one i € {1,2,..., n}. 

PROOF. We use induction on n. For the n = 1 base case the result is clear. 

For the inductive step, assume the inductive hypothesis: that the lemma 
holds for n such that 1 < n < k. We must show that it holds for n = k + 1. 
Assume that p is prime and that p \ a\a 2 ■ ■ ■ a^ak+i- Write 0102 • • • a*, as a, and 
dfc+i as b. Then p j a or p \ b by Lemma 7.3. If p \ a = a\ - ■ ■ then by the 
induction hypothesis, p \ ai for some i € {1,..., fc}. If p \ b then p \ ak+i- So 
we can say that p \ ai for some i € {1,2,... ,k + 1}. This verifies the lemma for 
n = k + 1. Hence by mathematical induction, it holds for all n > 1. QED 
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CHAPTER 7. THE FUNDAMENTAL THEOREM 


7.5 Lemma (Fundamental Theorem, Existence) If n > 1 then there exist 
primes pi, ■ ■ ■ ,p s , where s > 1, such that n = P 1 P 2 • • ■ p s and pi < p 2 < ■ ■ ■ < p s - 

PROOF. We will use induction on n. The base step is n = 2: in this case, since 
2 is prime we can take s = 1 and pi = 2. 

For the inductive step, assume the hypothesis that the lemma holds for 
2 < n < k: we will show that it holds for n = k + 1. If k + 1 is prime then s = 1 
and pi = k + 1. If k + 1 is composite then write k + 1 = ab where 1 < a < k + 1 
and 1 < b < k + 1. By the induction hypothesis there are primes pi,... ,p u and 
qi,... ,q v such that a = p\ ■ ■ -p u and b = q± - ■ ■ q v . This gives that k + 1 is a 
product of primes 

k + 1 = ab = pip 2 ■ ■ ■ p u qiq -2 ■■■q v , 

where s = u + v. Reorder the primes into ascending order, if necessary. 

The base step and the inductive step together give us that the statement is 
true for all n > 1. QED 

7.6 Lemma (Fundamental Theorem, Uniqueness) If n = P 1 P 2 • • -p s for 

s > 1 with pi < P 2 < • • • < p s , and also n = qiq 2 ■ ■ ■ qt for t > 1 with 
qi < 92 < • • • < qt, then t = s, and pi = qi for all i between 1 and s. 

PROOF. The proof is by induction on s. In the s = 1 base case, n = pi is prime 
and we have pi = qiq 2 • • • qt- Now, t must be 1 or else this is a factorization of 
the prime pi, and therefore pi = q-\. 

Now assume the inductive hypothesis that the result holds for all s with 
1 < s < k. We must show that the result then holds for s = k + 1. Assume that 
n = P1P2 • • • pkPk+i where pi < P2 < • ■ ■ < Pk+i, and also n = qiq2 ■ ■ ■ qt where 
qi < <?2 < • • • < qt- Clearly Pk+i \ n, so Pk+i \ qi • • • qt- Euclid’s Lemma then 
gives that Pk+i divides some qi. That implies that Pk+i = qt, or else Pk+i would 
be a non-1 divisor of the prime qt, which is impossible. Hence Pk+i = qi < qt- 
A similar argument shows that qt = Pj < Pk+i- Therefore Pk+i = qt- 
To finish, cancel Pk+i = qt. from the two sides of this equation. 

P 1 P 2 ■ ■ -PkPk+i = qiq 2 ■ ■ ■ qt-iqt 

Now the induction hypothesis applies: k = t — 1 and pi = qi for i = 1,... ,t— 1. 

So the lemma holds also in the s = k + 1 case, and so by mathematical 
induction it holds for all s > 1. QED 

7.7 Remark Unique factorization gives an alternative, conceptually simpler, 
way to find the greatest common divisor of two numbers. For example, 600 = 
2 3 • 3 1 • 5 2 • 7° and 252 = 2 2 • 3 2 • 5° • 7. Now, 2 3 divides both number. So does 
3 1 , but 3 2 does not divide both. Also, the highest power of 5 dividing both 
numbers is 5°, and similarly the highest power of 7 that works for both is 7°. 
So gcd(600, 252) = 2 2 • 3 1 • 5° • 7° = 24. In general, we can find the greatest 
common divisor of two numbers factoring, then taking the minimum power of 
2, times the minimum power of 3, etc. 
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The difficulty with this method is that we must factor the numbers. But 
factorization is very difficult! That is, for numbers that are large, factoring is 
slow while the Euclidean algorithm is relatively fast. 
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Chapter 8 

Distribution of Primes 


The Sieve of Eratosthenes is an ancient method to find primes. To find the 
primes less than n, list the numbers from 2 to n — 1. The smallest number, 2, is 
prime. Cross off all proper multiples of 2 (that is, the even numbers greater than 
2). The smallest number remaining, 3, is prime. Cross off all proper multiples of 
3, that is, 6, 9, etc. (some of them have already been eliminated). The smallest 
remaining number, 5, is prime. Cross off all proper multiples of 5. Continue 
this process until the list is exhausted. 

Here is what is left when the sieve filters out the nonprimes less than 100. 



00 01 

02 

03 04 

05 

06 07 08 

09 

0 


2 

3 

5 

7 


10 

11 


13 


17 

19 

20 



23 



29 

30 

31 




37 


40 

41 


43 


47 


50 



53 



59 

60 

61 




67 


70 

71 


73 



79 

80 



83 



89 

90 





97 



Obviously, the columns with even numbers and the columns with multiples of 
5 are empty (except for 2 and 5) but this is an artifact of the fact that the rows 
of the table are 10 = 2 • 5 wide. Other than that, at first glance no pattern is 
apparent. 


8.1 Theorem (Wilson’s Theorem) There are arbitrarily long gaps between 
primes: for any positive integer n there is a sequence of n consecutive composite 
integers. 

PROOF. Given n > 1, consider a = (n + 1)! + 2. We will show that all of the 
numbers a, a + 1,..., a + (n — 1) are composite. 
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Since n+ 1 > 2, clearly 2 | (n+1)!. Hence 2 | (n+l)!+2. Since (n+l)!+2 > 2, 
we therefore have that a = (n + 1)! + 2 is composite. We will finish by showing 
that the *-th number in the sequence, a + i where 0<z<n—1, is composite. 
Because 2<* + 2<n+l, we have that (i + 2) | (n + 1)!. Hence i + 2 | a + i = 
(n+l)! + (i+2). Because a+i > i+2 > 1, we have that a+i is composite, qed 

8.2 Definition For any positive real number x, the number of primes less than 
or equal to x is n(x). 

For example, 7r(10) = 4. 

The next result was first conjectured in 1793 by by Gauss, on the basis of 
numerical evidence like that in the table above. It was, however, not proved 
until over 100 years later, by Hadamard and Vallee Poussin. The proof is beyond 
the scope of this course. 

8.3 Theorem (The Prime Number Theorem) 


*^oo(a/ln(s)) 

Here is a table of values of 7r(10*) and 10*/ln(10 l ) for i = 2,..., 10 (the 
second set of values have been rounded to the nearest integer). 


X 

n(x) 

round (x/ ln(x)) 

10 2 

25 

22 

10 3 

168 

145 

10 4 

1229 

1086 

10 5 

9592 

8686 

10 6 

78498 

72382 

10 7 

664579 

620421 

10 8 

5761455 

5428681 

10 9 

50847534 

48254942 

10 10 

455052511 

434294482 


This table has been continued up to 10 21 , but mathematicians are still working 
on finding the value of 7r(10 22 ). Of course, computing the approximations are 
easy, but finding the exact value of 7r(10 22 ) is hard. 



Chapter 9 

Fermat Primes and 
Mersenne Primes 


A formula that produces the primes would be nice. Historically, lacking such 
a formula, mathematicians have looked for formulas that at least produce only 
primes. 

In 1640 Fermat noted that the numbers in this list 


n 

0 12 3 4 

F n = 2< 2 “> + 1 

3 5 17 257 65,537 


are all prime. He conjectured that F n is always prime. Numbers of the form 
2 2 +1 are called Fermat numbers. 

9.1 Lemma Let a > 1 and n > 1. If a n + 1 is prime then a is even and n = 2 k 
for some k > 1. 

PROOF. We first show that n is even. Suppose otherwise, and recall the well- 
known factorization. 

a n - 1 = (a - l)(a" -1 + a n ~ 2 + ■ ■ ■ + a + 1) 


Replace a by —a. 

(—a) n - 1 = (-o - 1) ((—a)"” 1 + (- a ) n ~ 2 + • • • + (-a) + l) 

If the exponent n is odd then n — 1 is even, n — 2 is odd, etc. So we have 
(— a) n = —a n , (—a) n_1 = a n ~ 1 , (—a)" -2 = —a" -2 , etc., and the factorization 
becomes 

-( a n + 1) = -{a + 1) (a™" 1 - a n ~ 2 +- a + l) . 

Then changing the sign of both sides gives ( a n + 1) = (a + l)(a” -1 — a n ~ 2 + 
• • • — a + 1). But with n > 2, we have l<a + l<a" + l. This shows that if n 
is odd and a > 1, then a n + 1 is not prime. 
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So n is even. Write n = 2 s ■ t where t is odd. Then if a n + 1 is prime we 
have (a 2 )* + 1 is prime. But by what we just showed this cannot be prime if t 
is odd and t > 2. So we must have t = 1 and therefore n = 2 s . 

Also, a n + 1 prime implies that a is even since if a is odd then so is a”, 
and in consequence a n + 1 would be even. But the only even prime is 2, adnd 
we are assuming that a > 1 and so we have a > 2, which implies that so 
a n + 1 > 3. QED 

9.2 Definition A prime number of the form F n = 2 < - 2 ^ + 1, n > 0, is a Fermat 
prime. 

Euler showed that Fermat number next on the table, F-j = 4, 294, 967,297, 
is composite. 

As n increases, the F n ’s increase in size very rapidly, and are not easy to 
check for primality. We know that F n is composite for all n such that 5 < n < 30, 
and a large number of other values of n including 382447 (the largest one that 
I know). Many researchers now conjecture that F n is composite for n > 5. So 
Fermat’s original thought that F n is always prime is badly mistaken. 

Mathematicians have also looked for formulas that produce many primes. 
That is, we can guess that numbers of various special forms are disproportion¬ 
ately prime. One form that has historically been of interest is are the Mersenne 
numbers M n — 2" — 1. 


n 

2 3 5 

7 

13 

f(n) 

3 7 31 

127 

8191 


All of the numbers on the second row are prime. Note that 2 4 — 1 is not prime, 
so this is not supposed to be a formula that gives only primes. 

9.3 Lemma Let a > 1 and n > 1. If a n — 1 is prime then a = 2 and n is prime. 

PROOF. Consider again a n — 1 = (a — l)(a n_1 + • • • + a + 1) Note that if a > 2 

and n > 1 then a — 1 > 1 and a 71-1 H-b a + 1 >a + l > 3 so both factors are 

greater then 1, and therefore a n — 1 is not prime. Hence if a" — 1 is prime then 
we must have a = 2. 

Now suppose 2” — 1 is prime. We claim that n is prime. For, if not, then 
n = st where 1 < s < n and 1 < t < n. Then 2™ — 1 = 2 st — 1 = (2 s ) 4 — 1 is 
prime. But we just showed that if a n — 1 is prime then we must have a = 2. So 
we must have 2 s = 2, and hence s = 1 and t = n. Therefore n is not composite, 
that is, n is prime. QED 

9.4 Corollary If M n is prime, then n is prime. 

PROOF. This is immediate from Lemma 9.3. QED 

At first it was thought that M p = 2 P — 1 is prime whenever p is prime. But 
in 1536, Hudalricus Regius showed that M\\ = 2 11 — 1 = 2047 is not prime: 
2047 = 23 • 89. 
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9.5 Definition A prime number of the form M n = 2 n — 1, n > 2, is a Mersenne 
prime. 

People continue to work on determining which Mp s are prime. To date 
(2003-Dec-09), we know that 2 p — 1 is prime if p is one of the following 40 
primes: 2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607, 1279, 2203, 2281, 
3217, 4253, 4423, 9689, 9941, 11213, 19937, 21701, 23209, 44497, 86243, 110503, 
132049, 216091, 756839, 859433, 1257787, 1398269, 2976221, 3021377, 6972593, 
13466917, and 20996011. 

The first number with more than a thousand digits known to be prime was 
Mi 253 - The largest number on that list was found on 2003-Nov-17. This number 
has 6, 320,430 digits. It was found as part of the Great Internet Mersenne Prime 
Search (GIMPS). (You can participate in this search by setting a program to 
run at times when your computer is not busy; see Chris Caldwell’s page for 
more about this.) Later we will see a connection between Mersenne primes and 
perfect numbers. 

One reason that we know so much about Mersenne primes is that the follow¬ 
ing test makes it easier to check whether or not M p is prime when p is a large 
prime. 

9.6 Theorem (The Lucas-Lehmer Mersenne Prime Test) Let p be an 

odd prime. Define the sequence r\ , rq, rq,. r p _i by the rules rq = 4, and for 
k> 2, 

rk = ( p fc_! - 2) mod M p . 

Then M p is prime if and only if r p _i = 0. 

The proof of this is beyond the scope of this book. 

9.7 Example Let p = 5. Then M p = M 5 =31. 

ri = 4 

7*2 = (4 2 — 2) mod 31 = 14 mod 31 = 14 
7*3 = (14 2 — 2) mod 31 = 194 mod 31 = 8 
7q = (8 2 — 2) mod 31 = 62 mod 31 = 0 

Hence by the Lucas-Lehmer test, M 5 = 31 is prime. 

9.8 Remark Note that the Lucas-Lehmer test for M p = 2 P — 1 takes only p— 1 
steps. On the other hand, if we try to prove that M p is prime by testing all 
primes less than or equal to \jM p then must consider about 2steps. This 
is much larger, in general, than p. 

No one knows whether there are infinitely many Mersenne primes. 
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Chapter 10 

The Functions a and r 


10.1 Definition Where n is a positive integer, r(n) is the number of positive 
divisors of n. 

10.2 Example The number 12 = 3 • 2 2 has positive divisors 1, 2,3,4, 6,12, and 
so r(12) = 6. 

10.3 Definition Where n is a positive integer, cr(n) is the sum of the positive 
divisors of n. 

A positive divisor d of n is a proper divisor if d < n. The sum of all proper 
divisors of n is <r*(n). 

Note that if n > 2 then cr*(n) = a(n) — n. 

10.4 Example cr(12) = 1 + 2 + 3 + 4 + 6 + 12 = 28, <r*(12) = 16. 

10.5 Definition A number n > 1 is perfect if a*(n ) = n. 

10.6 Example The first perfect number is 6 because its proper divisors are 1, 
2 and 3. 


10.7 Theorem Consider the prime factorization n = p^p^ 2 ■ ■ -Pr r ■ 


(1) r(n) 

(2) a(n) 


(ei + l)(e2 + 1) • • • ( e r + 1) 

^ 1 + 1 - 1 P? + 1 - 1 Pr r + 1 ~ 1 
Pi - 1 P'2 - 1 Pr - 1 


10.8 Example If n = 72 = 2 3 • 3 2 then r(72) = (3 + 1)(2 + 1) = 12 and 


*( 72 ) 


2 4 - 1 3 3 - 1 
2-1' 3-1 


15 • 13 = 195. 


Proof of item (1). From the Fundamental Theorem of Arithmetic, if d is a factor 
of n then the prime factors of d come from those of n. Hence d \ n iff d = 
Pi^pi 2 ''' Pr r where for each i, 0 < /* < e*. There are (ei + l)(e 2 +1) • • • (e r +1) 
choices for the exponents fi, / 2 ,..., f r - QED 
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Our proof of the second item requires two preliminary results. 

10.9 Lemma Suppose that n = ab, where a > 0, b > 0 and gcd(a, b) = 1. 
Then a(n) = a(a)a(b). 

PROOF. Since a and b have only 1 as a common factor, the Fundamental The¬ 
orem of Arithmetic, shows that d \ n only when d factors into d = d\d 2 where 
di | a and d 2 \ b. That is, the divisors of ab are products of the divisors of a 
with the divisors of b. Let the divisors of a be 1, ai,..., a s and the divisors of b 
be 1, b \,..., b t . These are the divisors of n = ab. 

1) bi, b 2 ,..., b t 

ai ■ l,ai • &i,ai • b 2 , ■ ■ ■,ai • b t 

a 2 ■ 1, a 2 ■ bi, a 2 • b 2 ,.. ■ , a 2 • b t 


a s ■ 1, a s ■ &i, a s ■ b 2 ,,.., a s ■ b t 

This list has no repetitions because, as gcd(a, b) = 1, if aibj = akbe then a. t = a^ 
and bj = be. Therefore to find cr(&) we can sum the rows 

1 + b\ H-+ b t = a(b) 

a il + ai&i + • • • + a\bt = aicr(6) 


a s • 1 + a s bi + • • • + a s b t = a s a{b) 


and add those partial sums 

a(n) = a(b) + aia(fo) + a 2 a(b ) + • • • + a 3 a(b) 
= (1 + di + a 2 + • • • + a s )cr(b ) 

= a(a)a(b) 

to get the required result. 

10.10 Lemma If p is a prime and k > 0 then 


a{p k ) 


p k+1 - 1 
p-1 


PROOF. Since p is prime, the divisors of p k are 1 ,p,p 2 ,.. ■ ,p k . Hence 

p k+i _ i 


a(p k ) = 1 +p + p 2 H- \-p k = 


p-1 


QED 


follows from the formula for the sum of a geometric series. 


QED 
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Proof of item (2). Let n = p^p^ 2 •••p® r . This proof is by induction on the 
number of prime factors r. In the r = 1 base case we have n = pf 1 and the 
result follows from Lemma 10.10. 

For the inductive step, the inductive hypothesis is that the statment is true 
when 1 < r < k. Consider the r = k + 1 case: n = p® 1 • • ■p t jfp e k k ff where the 
primes are distinct. Let a = p^ 1 ■ ■ ■ p k k and b = pjfff ■ Clearly gcd(a, b) = 1. 
Lemma 10.9 applies to give that a(n) = cr(a)a(b). The inductive hypothesis 
and Lemma 10.10 give 


cr(a) 


ei+1 -i 

Pi -i 
pi -1 


pl k+l -A 

Pk -1 J 


a(b) = 


pT + \ 1+1 i 

Pk +1 -1 


and therefore 


r(n) = 


pT +1 ~ 1 

Pi ~ 1 


( pl+i 1+l -A 

y Pk+i -1 J 


as desired. So the result holds for r = k + 1, and that implies that the theorem 
is true for all integers by the principle of mathematical induction. QED 



28 


CHAPTER 10. THE FUNCTIONS a AND r 



Chapter 11 

Perfect Numbers and 
Mersenne Primes 


A search for perfect numbers up to 10,000 finds only these. 

6 = 2-3 
28 = 2 2 • 7 
496 = 2 4 • 31 
8128 = 2 6 • 127 

Note that 3 = 2 2 — 1, 7 = 2 3 — 1, 31 = 2 5 — 1, and 127 = 2' — 1 are Mersenne 
primes. We can conjecture that all perfect numbers fit this pattern. This chapter 
discusses to what extent this is known to be true. 

11.1 Theorem If 2 P — 1 is a Mersenne prime then 2 P_1 • ( 2 P — 1) is perfect. 

PROOF. Write q = 2 P —1 and n = 2 p ~ 1 q. Since q is odd and prime, Theorem 10.7 
gives that <r(n) is 

u (2 p-1 g) = = (2 p - 1)(« + 1) = (2 P - l)2 p = 2n. 

That is, er(n) = 2 n, and so n is perfect. QED 

11.2 Theorem If n is even and perfect then there is a Mersenne prime 2 p — 1 
such that n = 2 P ” 1 (2 P — 1). 

PROOF. Suppose that n is even and perfect. Factor out all of the 2’s to get 
n = 2 k ■ q with q an odd number, and k > 1 since n is even. Since q is odd, 
gcd(2 fe ,g) = 1 and so by Lemmas 10.9 and 10.10 we have a(n) = a(2 k )a(q) = 
(2 fc+1 — l)cr(q). Thus, as n is perfect, 

2 k+1 q = 2 n = a(n) = (2 fc+1 - l)a(q). 
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Now substituting a(q) = a*(q) + q , into the prior displayed equation gives 
2 fc + 1 g=(2 fe+1 ^l )(<J*(q) + q) 

that is 

2 k+1 q = (2 k+1 - 1 )<r*(q) + 2 k+1 q - q 

This implies that 

(*) a*(q)(2 k + 1 -l) = q. 

So <J*(q) is a divisor of q. Since k > 1 we have that 2 k+1 — 1 > 4 — 1 = 3. So 
cr*(q) is a proper divisor of q. But cr*(q) is the sum of all of the proper divisors 
of q. This can only happen if q has only one proper divisor, that is, it implies 
that q is prime and so cr*(q) = 1. Then equation (*) shows that q = 2 fc+1 — 1. 
So q is a Mersenne prime and k + 1 = p is prime. Therefore n = 2 P ~ 1 • ( 2 P — 1), 
as desired. QED 

11.3 Corollary There is a one-to-one correspondence between the even perfect 
numbers and the Mersenne primes. 

Here are two questions that remain open: (i) Are there infinitely many even 
perfect numbers? (ii) Are there any odd perfect numbers? (We know that if an 
odd perfect number exists, then it must be greater than 10 50 .) 



Chapter 12 

Congruences 


12.1 Definition Let m > 0. We we say that the numbers a and b are congruent 
modulo to, denoted a = b (mod to), if a and b leave the same remainder when 
divided by to. The number in is the modulus of the congruence. The notation 
a jk b (mod to) means that they are not congruent. 

12.2 Lemma The numbers a and b are congruent modulo to if and only if 
to | (a — b), and also if and only if m\ (b — a). 

PROOF. Write a = mq a + r a and b = mqb + rb for some q a , qb, r ai and 77 , , with 
0 < r 0 , rh < to. Subtracting gives a — b = m(q a — qb) + (r 0 — rb). Observe that 
the restrictions on the remainders imply that —to < r a — n < to, and so r a — rb 
is not a multiple of to unless r a — rb = 0. 

If a and b are congruent modulo to then r a = rb, which implies that a — b = 
m(q a — qb), which in turn gives that a — b is a multiple of to. 

The implications in the prior paragraph reverse: if a — b is a multiple of m 
then in the equation a — b= m(q a — qb) + ( r a — rb) we must have that r a — H = 0 
by the observation in the first paragraph, and therefore r a = rb. 

The b — a statement is proved similarly. QED 

12.3 Examples 

1. 25 = 1 (mod 4) since 4 | 24 

2. 25 ^ 2 (mod 4) since 4 \ 23 

3. 1 = —3 (mod 4) since 4 | 4 

4. a = b (mod 1) for all a, b 

5 . a = b (mod 0) a = 6 for all a, b 

Do not confuse the use of mod in Definition 12.1 
a = b (mod m) if to | a — b 
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with that of Definition 3.3. 

a mod b = r where r is the remainder when a is divided by b 
The two are related but not identitical. 

12.4 Example One difference between the two is that 25 = 5 (mod 4) is true 
while 25 = 5 mod 4 is false (it asserts that 25 = 1). 

The ‘mod’ in a = b (mod m) defines a binary relation, a relationship between 
two things. The ‘mod’ in a mod b is a binary operation, just as addition or 
multiplication are binary operations. Thus, 

a = b (mod to) •£=> a mod to = b mod m. 

That is, if m > 0 and a = r (mod m ) where 0 < r < m then a mod m = r. 
Expressions such as 


x = 2 
4 2 = 16 

x 2 + 2x = sin(x) + 3 

are equations. By analogy, expressions such as 

x = 2 (mod 16) 

25 = 5 (mod 5) 
x 3 + 2x = Qx 2 + 3 (mod 27) 


are called congruences. 

The next two theorems show that congruences and equations share many 
properties. 

12.5 Theorem Congruence is an equivalence relation: for all a, &, c, and 
to > 0 we have 

(1) (Reflexivity property) a = a (mod m) 

(2) (Symmetry property) a = b (mod to) => b= a (mod to) 

(3) (Transitivity property) a = b (mod to) and b = c (mod to) => a = c 
(mod to) 

PROOF. For reflexivity: on division by to, any number leaves the same remainder 
as itself. 

For symmetry, if a leaves the same remainder as 6, then b leaves the same 
remainder as a. 

For transitivity, assume that a leaves the same remainder as b on division 
by to, and that b leaves the same remainder as c. The all three leave the 
same remainder as each other, and in particular a leaves the same remainder as 
c. QED 
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Below we will consider polynomials f(x) = a n x n +a n -\x n ~ 1 + ■ ■ ■ + aiX + ao- 
We will assume that the coefficients a n ,...,a o are integers and that x also 
represents an integer variable. Here the degree of the polynomial is an integer 
n > 0 . 

12.6 Theorem If a = b (mod to) and c = d (mod to), then 

( 1 ) a + c = b + d (mod to) and a — c = b — d (mod to) 

( 2 ) ac = bd (mod to) 

(3) a n = b n (mod m) for all n > 1 

(4) /(a) = f(b) (mod m ) for all polynomials f(x) with integer coefficients. 

Proof of (1). Since a — c = a + (—c), it suffices to prove only the addition case. 
By assumption m \ a — b and to | c — d. By linearity of the ‘divides’ relation, 
to | (a — b) + (c — d), that is in \ (a + c) — (b + d). Hence a + c = b + d 
(mod m). QED 

Proof of (2). Since m \ a — b and to | c — d, by linearity m \ c(a — b) + b(c — d). 
Now, c(a — b) + b(c — d) = ca— bd , hence m \ ca— bd , and so ca = bd (mod to), 
as desired. QED 

Proof of (3). We prove this by induction on n. If n = 1, the result is true by 
the assumption that a = b (mod m ). Assume that the result holds for n = 
1,... ,k. Then we have a k = b k (mod to). This, together with a = b (mod m ) 
using property (2) above, gives that aa k = bb k (mod to). Hence a k+1 = b k+1 
(mod m) and the result holds in the n = k + 1 case. So the result holds for all 
n > 1, by induction. QED 

Proof of (f). Let f(x) = c n x n + ■ ■ ■ + C\X + Cq. We prove by induction on 
the degree of the polynomial n that if a = b (mod to) then c n a n + ■ ■ ■ + Co = 
c n b n + ■■■ + cq (mod to). For the degree n = 0 base case, by the reflexivity of 
congruence we have that Co = Co (mod to). 

For the induction assume that the result holds for n = k. Then we have 

(*) Cka k + • • • + cia + Co = Ckb k + • • • + C\b + Cq (mod m ). 

By item (3) above we have a k+1 = b k+1 (mod to). Since Ck+\ = c^+i (mod to), 
using item ( 2 ) above we have 

(**) c k+1 a k+1 = c k+1 b k+1 (mod to). 

Now we can apply Theorem 15.3 (1) to (*) and (**) to obtain 

Ck+ia k+1 + c k a k + -1- c 0 = c k+ \b k+1 + c k b k + -b c 0 (mod to). 

So by induction the result holds for all n > 0. QED 

12.7 Example (From [1].) The first five Fermat numbers 3, 5, 17, 257, and 
65, 537 are prime. We will use congruences to show that F 5 = 2 32 +1 is divisible 
by 641 and is therefore not prime. 
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Everyone knows that 2 2 = 4, 2 4 = 16, and 2 8 = 256. Also, 2 16 = (2 8 ) 2 = 
256 2 = 65,536. A straightforward division shows that 65,536 = 154 (mod 641). 

Next, for 2 32 , we have that (2 16 ) 2 = (154) 2 (mod 641). That is, 2 32 = 
23,716 (mod 641). Since an easy division finds that 23,716 = 640 (mod 641), 
and 640 = —1 (mod 641), we have that 2 32 = —1 (mod 641). Hence 2 32 + l = 0 
(mod 641), and so . 641 | 2 32 + 1, as claimed. Clearly 2 32 + 1 yf 641, so 2 32 + 1 
is composite. 

The work done here did not require us to find the value of 2 32 + 1 = 
4,294,967,297 and divide it by 641; instead the calculations were with much 
smaller numbers. 



Chapter 13 

Divisibility Tests 


Elementary school children know how to tell if a number is even, or divisible by 
5, by looking at the least significant digit. 

13.1 Theorem If a number a has the decimal representation a„_ 1 IO " -1 + 
a„_ 2 10” -2 + • • • + ailO + ao then 

(1) a mod 2 = ao mod 2 

(2) a mod 5 = ao mod 5 

PROOF. Consider f(x) = a n _ix™” 1 + - • - + ai;r + ao. Note that 10 = 0 (mod 2). 
So by Theorem 12.6 

an-ilO" T • • • T ailO T ao = a n _i0 n aiO T ao (mod 2). 

That is, a = a 0 (mod 2); this proves item (1). Since 10 = 0 (mod 5) also, the 
proof of item (2) is similar. QED 

13.2 Example Thus, the number 1457 is odd because 7 is odd: 1457 mod 2 = 
7 mod 2 = 1. And on division by 5 it leaves a remainder of 1457 mod 5 = 
7 mod 5 = 2. 

13.3 Theorem Where a — On-ilO”^ 1 + a n _ 2 10" -2 + • • • + ailO + a 0 is the 
decimal representation, 

(1) a mod 3 = ( a n -i + • • • + ao) mod 3 

(2) a mod 9 = (a„_i + • • • + ao) mod 9 

(3) a mod 11 = (ao — a\ + a 2 — + ■ ■ ■) mod 11. 

PROOF. Note that 10 = 1 (mod 3). Theorem 12.6 gives 

a„_ilO n • • • “t - ailO ao = a n _il n 4- • • • -t- ail -t- ao (mod 3) 

and so a = a n -\ + ■ ■ ■ + ai + ao (mod 3). This proves item (1). Since 10 = 1 
(mod 9) also, the proof of item (2) is similar. 
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For item (3), note that 10 = —1 (mod 11) so 
a„_il0" 1 + • • • + cti 10 + ciq = a n _i (—1)” 1 + • • • + ai(—1) + ag (mod 11). 
That is, a = ao — ai + a 2 — • • • (mod 11). QED 

13.4 Example Consider 1457 again. For divisibility by 3 we have 1457 mod 

3 = (1 + 4 + 5 +7) mod 3 = 17 mod 3 = 8 mod 3 = 2. As for 9, we get 

1457 mod 9 = (1 + 4 + 5 +7) mod 9 = 17 mod 9 = 8 mod 9 = 8. Finally, for 11, 

the calculation is 1457 mod 11 = 7 — 5 + 4 — 1 mod 11 = 5 mod 11 = 5. 

Note that m \ a <+> a mod m = 0 so from the prior two results we obtain 
immediately the following. 

13.5 Corollary Let a = a n _ilO” -1 + a„_ 2 10 n_2 + • • • + ailO + ao- 

(1) 2 | a <+> ao = 0,2,4, 6 or 8 

(2) 5 | a <+> ao = 0 or 5 

(3) 3 | a +> 3 | ao + ai + • • • + a n _ l 

(4) 9 | a •+> 9 | ao + ai + • • • + a n _ i 

(5) 11 | a <+> 11 | a 0 — ai + a 2 — a 3 + • • •. 

13.6 Theorem Let a = a r !0 r + • • • + a 2 10 2 + ailO + ao be the decimal repre¬ 
sentation, so that we write a as the sequence a r a r -i ■ • • aiao- Then 

(1) 7 | a <+> 7 | a r ■ ■ ■ a\ — 2ao- 

(2) 13 | a <+> 13 | a r • • • ai — 9ao 

(where a r • • • ai is the sequence representing (a — ao)/10). 

PROOF. For item (1), let c = a r ■ ■ ■ ai so that a = 10c+ao. Since gcd(7, —2) = 1 
we have that 7 | a <+> 7 | —2a. Consequently, consider —2a = —20c — 2ao- 
Because 1 = —20 (mod 7), we have that —2a = c — 2ao (mod 7). Therefore, 
7 | —2a <t=> 7 | c — 2ao- It follows that 7 | a <t=> 7 | c — 2ao, which is what we 
wanted to prove. 

The proof of item (2) is similar. QED 

13.7 Example We can test whether 7 divides 2481. 

7 | 2481 <+> 7 | 248 - 2 <+> 7 | 246 <+ 7 | 24 - 12 <+ 7 | 12 
Since 7 j 12 we have that 7 j 2481. 

13.8 Example The number 12987 is divisible by 13 because 

13 | 12987 +> 13 | 1298 - 63 +> 13 | 1235 +> 13 | 123 - 45 +> 13 | 78 


and 13 • 6 = 78. 




Chapter 14 

More Properties of 
Congruences 


Theorem 12.6 provides some laws of algebra for =. A typical algebra problem is 
to solve for an unknown; for instance, we can look for x such that 2x = 7 mod 15. 

14.1 Theorem Let m > 2. If a and m are relatively prime then there exists 
a unique integer a* such that aa* = 1 (mod m) and 0 < a* < m. 

PROOF. Assume that gcd(a, m) = 1. Bezout’s Lemma applies to give an s and 
t such that as + mt = 1. Hence as — 1 = that is, m | as — 1 and so 

as = 1 (mod m). Accordingly, let a* = s mod m so that 0 < a* < m. Then 
a* = s (mod m) so aa* = 1 (mod m). 

To show uniqueness, assume that ac = 1 (mod in ) and 0 < c < m. Then 
ac = aa* (mod m). Multiply both sides of this congruence on the left by c and 
use the fact that ca = 1 (mod m ) to obtain c = a* (mod m). Because both are 
in [0 .. m), it follows that c = a*. QED 

We call a* the inverse of a modulo m. Note that we do not denote a* by 
a -1 here since we keep that symbol for the usual meaning of inverse. 

14.2 Remark The proof shows that Blankinship’s Method will compute the 
inverse of a, when it exists. But for small m we may find a* by trial and error. 
For example, take m = 15 and a = 2. We can check each possibility: 2-0^1 
(mod 15), 2-1^1 (mod 15), ..., 2-8 = 1 (mod 15). So we can take 2* = 8. 

Note that we may well have ca = 1 mod m with c yf a if c = a* (mod m) 
and c > m or c < 0. For instance, 8-2=1 mod 15 and also 23 • 2 = 1 mod 15. 
So the inverse is unique only if we specify that 0 < a* < m. 

The converse of Theorem 14.1 holds. 

14.3 Theorem Let m >0. If ab = 1 (mod m ) then both a and b are relatively 
prime to m. 
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PROOF. If ab = 1 (mod to), then to \ ab — 1. So ab — 1 = mt for some t. Hence, 
ab + m(—t) = 1. 

The proof of Bezout’s Lemma, Lemma 5.3, shows that gcd(a, to) is the small¬ 
est positive linear combination of a and to. The last paragraph shows that 
there is a combination that adds to 1. Since no combination can be positive 
and smaller than 1, we have that gcd(a,m) = 1. The case of gcd(6, to) is 
similar.. QED 

14.4 Corollary A number a has an inverse modulo to if and only if a and to 
are relatively prime. 

The second paragraph of Theorem 14.1 uses a technique that is worth iso¬ 
lating. 

14.5 Theorem (Cancellation) Let to > 0. If gcd(c, to) = 1 then ca = cb 
(mod to) => a = b (mod to). 

PROOF. If gcd(c, to) = 1 then it has an inverse c* modulo to, such that c*c = 1 
(mod to). Since ca = cb (mod to) by Theorem 12.6, c*ca = c*cb (mod to). But 
c*c = 1 (mod to) so c*ca = a (mod to) and c*cb = b (mod to). By reflexivity 
and transitivity this yields a = b (mod m). QED 

Although in general we cannot cancel if gcd(c, m) > 1, the next result is 
some consolation. 

14.6 Theorem If c > 0 and to > 0 then a = b (mod to) -o- ca = cb (mod cm). 

PROOF. The congruence a = b (mod to) is true if and only if to | (a — b) holds, 
which in turn holds if and only if cm \ (ca — cb). QED 

14.7 Theorem Fix to > 0 and let d = gcd(c, to). Then ca = cb (mod to) => 
a = b (mod m/d). 

PROOF. Since d = gcd(c, to), the equations c = d(c/d) and to = d(m/d) involve 
integers. Rewriting ca = cb (mod to) gives 

d (/) asd (/) b ( m ° dd (-j)y 

By Theorem 14.6 we have 

G) os ()) i ' (mod 7>- 

Since d = gcd (c, to), we have that gcd (c/d,m/d) = 1 and so by cancellation, 
Theorem 14.5, a = b (mod m/d). QED 

14.8 Theorem If to > 0 and a = b (mod to) then gcd(a,m) = gcd(6, to). 

PROOF. Let d a = gcd (to, a) and = gcd (to, 6). Since a = b (mod to) we have 
a — b = mt for some t. Rewrite that as a = mt + b and note that db \ m and 
db | b , so db | a. Thus, db is a common divisor of to. and a, and so db < d a ■ A 
similar argument gives that d a < db, and therefore (4 = d a - QED 
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14.9 Corollary Fix m > 0. If a = b (mod m ) then a has an inverse modulo 
m if and only if b does also. 

PROOF. Immediate. QED 
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Chapter 15 

Residue Classes 


The work that we’ve seen shows that if a = b (mod to) then the two numbers 
a and b, while not necessarily equal, are in some ways alike. 

15.1 Definition Fix m > 0. The residue class class of a modulo to (or 
congruence class, or equivalence class of a modulo m) is [a] = {x \ x = a 
(mod m)}, the set of all integers congruent to a modulo to. 

Note that, by definition, [a] is a set. 

[a] = {mq + a \ q £ Z} = {..., —2 in + a, —m + a, a, m + a, 2 to + a,... } 

Note also that [a] depends on to and so it would be more accurate to write [a\ m 
instead, but this would be cumbersome. 

15.2 Theorem If to > 0 then [a] = [6] <t=> a = b (mod m). 

PROOF. First assume that [a] = [6], Note that a £ [a] because a = a (mod to). 
And, because [a] = [6], we have a £ [ b ]. By definition of [6], then a = b (mod to). 

For the implication the other way, assume that a = b (mod to), aiming to 
prove that the sets [a] and [6] are equal. To prove that the sets are equal, we 
will prove that every element of the first is a member of the second, and vice 
versa. Suppose that x £ [a], so that x = a (mod to). Since a = b (mod to), by 
transitivity of equivalence, x = b (mod to), and so x € [6]. The argument to 
show that ii x £ [6] then x £ [a] is similar. QED 

15.3 Theorem Given to > 0. For every a there is a unique r £ [0..?n) such 
that [a] = [r]. 

PROOF. Let r = a mod in so that 0 < r < to, and a = r (mod to), and by 
Theorem 15.2, [a] = [r]. To prove that r is unique, suppose that [a] = [?’'], where 
0 < r' < to. By Theorem 15.2, this implies that a = r' (mod to). This, together 
with the restriction that 0 < r' < to, implies that r' — a mod to = r. QED 

15.4 Theorem Given to > 0, there are exactly to distinct residue classes 
modulo to, namely, [0], [1],..., and [to — 1]. 


41 



42 


CHAPTER 15. RESIDUE CLASSES 


PROOF. By Theorem 15.3 we know that every residue class [a] is equal to one 
of [0], or [1], ..., or [to — 1]. So any residue classes is in this list. These residue 
classes are distinct: if 0 < rq < m and 0 < < m and [n] = \r^ then by the 

uniqueness part of Theorem 15.3 we must have r± = r%. QED 

15.5 Definition Any element x € [a] is a class representative. The element of 
[a] that is in [0.. to) is the principle class representative or principle residue. 



Chapter 16 

Z m and Complete Residue 
Systems 


Throughout this section we assume a fixed modulus m > 0. 

16.1 Definition The set {[a] \ a G Z} of all residue classes modulo m is 
denoted Z m . 

Recall that in a set, the order in which elements appear does not matter, 
and repeat elements collapse: the set { 0 , 2 ,3,1} and the set { 2 , 0 , 2 ,3,1,4,1} are 
equal. So, while at first glance Z m may seem to have infinitely many elements 
Z TO = {..., [—2], [—1], [0], [1], [2],... }, Theorem 15.4 shows that after the repeats 
collapse Z m = {[0], [1],..., [m — 1]}, and so instead Z m has exactly m elements. 

16.2 Example Fix m = 4. Then [1] = {..., —7, —3,1,5,...}, and so all of 
these classes are equal: • • • = [—7] = [—3] = [1] = [5] = • • •. We could, therefore, 
instead of Z 4 = {[0], [1], [2], [3]}, write Z 4 = { [ 8 ], [5], [— 6 ], [11]}. 

16.3 Definition A set of m integers {ao, cq,..., a m _i} is a complete residue 
system modulo m (or a complete set of representatives for Z m ) if the set Z m 
equals the set {[a 0 ], [ax ],..., [a m _ 1 ]}. 

16.4 Example These are complete residue systems modulo 5. 

1. {0,1, 2, 3,4} 

2 . {- 2 ,- 1 , 0 , 1 , 2 } 

3. {-9,14,12,10,8} 

4. {0 + 5ni, 1 + 5 ri 2 , 2 + 5713 ,3 + 5 ri 4 ,4 + 5774 }, where rq, 712 , 713 , 714,715 may 
be any integers. 

For each m > 0 there are infinitely many distinct complete residue systems 
modulo in. In particular, {0,1,..., m— 1} is the set of least nonnegative residues 
modulo m. 
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16.5 Theorem Fix m > 0. If to = 2k then {0,1,2,...,/:— 1, k,—(k — 
1), ..., —2, —1} is a complete residue system modulo m. If m = 2k + 1, then 
{0, 1,2 ,k, —k, . .., —2, —1} is a complete residue system modulo to . 

PROOF. If to = 2k, then since Z m = {[0], [1],..., [fe], [fc+l],..., [k+i], [k+k— 1]}, 
it suffices to note that [k + z] = [k + i — 2k] = {—k + i] = [— (k — *)]. So 
[k + 1] = [— (k — 1)], [k + 2] = [— (k — 2)],..., [k + k — 1] = [—1], as desired. 

In the n = 2fc+l case, [k+i] = [— (2k+l) + k+i] = [— k+i+1] = [— (k— i+1)] 
so [k + 1] = [—k], [k + 2] = [— (k — 1)],..., [2k] = [—1], as desired. QED 

16.6 Definition The complete residue system modulo in given in the prior 
theorem is the least absolute residue system modulo to. 

16.7 Example Where m = 2 32 , the least absolute residue system is 

{—(2 31 - 1), —(2 31 - 2),..., -2, -1,0,1, 2,..., 2 31 }. 



Chapter 17 

Addition and Multiplication 
in Z m 


In this chapter we show how to define addition and multiplication of residue 
classes modulo m. With respect to these binary operations Z m is a ring as 
defined in Appendix A. 

17.1 Definition For [a], [6] £ Z m , the sum of the residue class [a] and the 
residue class [6] is the residue class [o + b]. The product of the residue class [a] 
and the residue class [6] is the residue class [ab\. That is, 

[a] + [6] = [a + b] [a] [6] = [ab]. 

17.2 Example For m = 5 we have [2] + [3] = [5] and [2] [3] = [6]. Note that 
since 5 = 0 (mod 5) and 6 = 1 (mod 5) we can also write [2] + [3] = [0] and 

[2][3] = [1]. 

We must check that these binary operations are well defined. That is, since 
a residue class can have many representatives, we must check that the results 
of an operation do not depend on the representatives chosen for that operation. 

For example, fix m = 5 and consider [7] + [11]. We know that the residue 
classes [7] and the residue class [2] are equal, and also that [11] = [21]. Therefore 
for the binary operations to make sense we must have that [7] + [11] = [2] + [21]. 
In this case, [7] + [11] = [18] and [2] + [21] = [23], and [18] = [23] so this one 
example is fine. 

17.3 Theorem The results of the sum and product of residue classes does not 
depend on the choice of class representative: for any modulus m > 0, if [a] = [&] 
and [c] = [d] then [a] + [c] = [ b} + [d] and [a] [c] = [b] [d]. 

PROOF. This follows immediately from Theorem 12.6. QED 

When performing addition and multiplication in Z TO , we may at any time 
change class representatives, rewriting [a] by [a'], where a = a' (mod m). 
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17.4 Example Take m = 151 and consider the calculation [150] [149]. Then 
150 = —1 (mod 151) and 149 = —2 (mod 151), and so [150][149] = [—1][—2] = 
[ 2 ], an easier calculation. 


When working with Z m it is often useful to write all residue classes in the 
least nonnegative residue system, as we do in constructing the following addition 
and multiplication tables for Z 4 . 


[ 0 ] 

[ 1 ] 

[ 2 ] 

[3] 


[0] [1] [2] [3] 

[0] [1] [2] [3] 

[1] [2] [3] [0] 

[2] [3] [0] [1] 

[3] [0] [1] [2] 


[ 0 ] 

[ 1 ] 

[ 2 ] 

[3] 


[0] [1] [2] [3] 

[ 0 ] [ 0 ] [ 0 ] [ 0 ] 

[0] [1] [2] [3] 

[ 0 ] [ 2 ] [ 0 ] [ 2 ] 

[0] [3] [2] [1] 


Notice that we have reduced results of the sum and product to keep the repre¬ 
sentative in [0 .. 4). That is, in constructing those tables we follow the alogrithm 
that resclassa + [ 6 ] = [(a + b) mod m] and [a] [5] = [(a&) mod m\. 

This leads to an alternative way to define Z m and addition and multiplication 
in Z m . For clarity we will use different notation. 

17.5 Definition For m > 0, let J m be the set = {0,1,2,.,., m — 1} endowed 
with two binary operations: for a, b £ J m , let a © b = (a + b) mod m and 
a © b = {ab) mod m. 

Here are the addition and multiplication tables for J 4 . 


© 

0 

1 

2 

3 


0 12 3 

0 12 3 

12 3 0 

2 3 0 1 

3 0 12 


© 

0 

1 

2 

3 


0 12 3 
0 0 0 0 
0 12 3 
0 2 0 2 
0 3 2 1 


17.6 Remark The precise expression of the intuition that J m with © and 
0 is just like Z m with addition and multiplication is to say that the two are 
“isomorphic.” In this book we will leave the idea as informal. 


17.7 Example Let’s solve the congruence 272x = 901 (mod 9). Using residue 
classes modulo 9 we see that this congruence is equivalent to [272a;] = [901], 
which is equivalent to [272][a;] = [901]. That is equivalent to [2][x] = [1]. We 
know [a;] € {[0], [1],..., [ 8 ]}, so by trial and error we see that x = 5 is a solution. 



Chapter 18 

The Group of Units 


This is the multiplication table for Zg. 


© 

[ 0 ] 

[ 1 ] 

[ 2 ] 

[3] 

[4] 

[5] 


[0] [1] [2] [3] [4] [5] 


[ 0 ] [ 0 ] [ 0 ] [ 0 ] [ 0 ] [ 0 ] 

[0] [1] [2] [3] [4] [5] 

[0] [2] [4] [0] [2] [4] 

[0] [3] [0] [3] [0] [3] 

[0] [4] [2] [0] [4] [2] 

[0] [5] [4] [3] [2] [1] 


Note that some rows, and some columns, contain all of the members of Zg 
while others do not. We can state that as: for some [a], [ 6 ] £ Z 6 the equation 
[a] © x = [ 6 ] has no solution x. 

18.1 Example The equation [5] Ox = [3] has the solution x = [3]. In fact, for 
any [6] £ Z 6 , the equation [5] © x = [b] has a solution. However, the equation 
[4] © x = [1] has no solution. 

18.2 Definition Let m > 0. A residue class [a] £ Z TO is a unit if there is 
another residue class [b] £ Z TO such that [a] © [ b} = [1], In this case [a] and [b] 
are said to be inverses of each other in Z m . 

18.3 Theorem Let m > 0. A residue class [a] £ Z m is a unit if and only if 
gcd(a, m) = 1 . 

PROOF. Let [a] be a unit. Then there is a [b] such that [a] © [ 6 ] = [1]. Hence 
[ah] = [1] and so ab = 1 (mod m). Thus, by Theorem 14.3, gcd(a,?n) = 1. 

To prove the converse, let gcd(a, m) = 1. By Theorem 14.1 there is an 
integer a* such that aa* = 1 (mod m). Hence [aa*] = [1]. So [a] © [a*] = [1] 
and we can take b = a*. QED 


Note that from Theorem 14.8 if [a] = [ 6 ] -that is, if a = b (mod m) —then 
gcd(a,m) = 1 <t=> gcd( 6 , m) = 1. So, in checking whether or not a residue class 
is a unit we can use any representative of the class. 
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18.4 Theorem For m > 0, the set of units in Z m is the set of residue classes 
{[*] | 1 < * < to and gcd(*,?n) = 1 }. 

PROOF. If [a] € Z m then [a] = [i], where 0 < i < m — 1, so for each m > 0 we 
need only consider residue classes with representatives in the interval [0 ..to). 

If to = 1 then Z m consists of a single residue class Zi = {[0]} = {[1]}. Since 
[ 1 ] © [ 1 ] = [ 1 ], we have that this single class [ 1 ] is a unit. 

If to > 1 then gcd(0, to) = m ^ 1 and gcd(m, to) = to ^ 1, but gcd(i, to) = 1 
for 1 < i < m. So the theorem follows from Theorem 18.3. QED 

18.5 Definition The set of all units in Z m , the group of units, is denoted U m . 
(See Appendix A for the definition of a group.) 

18.6 Example Here are the first few U m ’s. 


i 

1 2 3 4 5 6 

Ui 

{[!]} {Ml {[1],[2]} {[1], [3]} {[1],[2],[3],[4]} {[1], [5]} 


18.7 Theorem The set of units U rn has these properties. 

1. (Closure) If [a] and [b] are members of U m then the product [a] [ 6 ] is also 
a member of U m . 

2 . (Associativity) For all [a], [ 6 ], [c] in U m we have ([a] © [ 6 ]) © [c] = [a] 0 

3. (Existence of an identity) [1] © [a] = [a] © [1] = [a] for all [a] € U m . 

4. (Existence of inverses) For each [a] £ U m there is a [a]* € U m such that 

M © H* = I 1 ]- 

5. (Commutativity) For all [a], [b] £ U m , we have that [a] © [b] = [b] © [a], 

18.8 Example Theorem 18.3 shows that 

U 15 = {[1], [2], [4], [7], [ 8 ], [11], [13], [14]} 

= {[!]. [2], [4], [7], [-7], [-4], [-2], [-1]}. 

Rather than list the entire multiplication table, we just show the inverse of each 
element. 


r 

[ 1 ] 

[ 2 ] 

[4] 

[7] 

[ 8 ] 

[ 11 ] 

[13] 

[14] 

inverse of r 

[ 1 ] 

[ 8 ] 

[4] 

[13] 

[ 2 ] 

[ 11 ] 

[7] 

[14] 


18.9 Theorem Let to > 0 and fix [a], [ 6 ] £ U m . Then the equation [a] 0 x 
resclassb has a unique solution x £ U m . 

PROOF. To see that it has a solution, consider [a]* ©[6]. By the closure property, 
that is an element of U m . Also, [a] 0 ([a]* © [b]) = ([a] 0 [a]*) © [ b } = [1] © [6] = 
[1 • b] = [b] , as required (the first equality follows by the associative property). 

To see that the solution is unique, suppose that x, x' £ U m are such that 
[a] © x — [ 6 ] and also [a] © x' = [ b ]. Then [a] © x = [a] © x'. Multiplying both 
sides of that equation by the inverse [a]* gives [a]* © ([a] ©x) = ([a]* © [a]) Qx = 
[1] © x = x on the left, and x' on the right. So the two are equal. QED 
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18.10 Definition If X is a set, the cardinality |X| is the number of elements 
in X. 

18.11 Example |{1}| = 1, |{0,1, 3,9}| = 4, |Z m | = m if m > 0. 

18.12 Definition If to > 1 then the Euler phi function (or the totient ) is 
0(m) — |{i € Z | 1 < i < m and gcd(z, m) = 1}|. 

18.13 Example Here are the first few values of 0. 


i 

1 

2 3 4 5 6 


1 

1 2 2 4 2 


Compare this to the table in Example 18.6. 

18.14 Corollary If m > 0 then U rn = (j>(m). 

Note that if p is any prime then 0(p) = p — 1. 

In general, though, 0(m) is not easy to calculate. However, computing 0(m) 
is easy once we know the prime factorization of m. 

18.15 Theorem Fix a, 6 > 0. If gcd(a, b) = 1 then 0(a6) = 0(a)0(6). 

18.16 Theorem If p is prime and n > 0 then 0 ( p n ) = p n — p n_1 . 

18.17 Theorem Let p!,p 2 , ■ ■ ■ ,Pk be distinct primes and let rii,ri 2 ,... , be 
positive integers. Then 

0 (pTpT ■■■pT) = ( pT -p?- 1 ) ■ ■ ■ (pT - pT~ 1 ) ■ 

The proofs of Theorem 18.15 and Theorem 18.17 are routine arguments by 
induction on n, and are left as exercises. 

Proof of Theorem 18.16. We want to count the number of elements in the set 
A = {1,2 ,,p n } that are relatively prime to p n . Let B be the set of elements 
of A that are not relatively prime, that is, that have a factor greater than 1 in 
common with p n . The nuber p is prime, so the only factors of p n are 1 ,p,... ,p n , 
and hence b = pk for some k. It follows that if a number b is an element of B 
then it has the form b = kp for some 1 < k < p™ -1 . That is, B is a subset of 
this set: {p, 2 p, 3 p,...,kp,... ,p n ~ 1 p}. But obviously every element of that set 
is not relatively prime to p n , so in fact B equals that set. 

The number of elements in A is H = p n and the number in B is \B\ = p" -1 , 
so the number of elements of A that are not in B is p n — p" -1 . QED 

18.18 Example 0(12) = 0(2 2 • 3) = (2 2 - 2 1 )(3 1 - 3°) = 2 • 2 = 4 

18.19 Example 0(9000) = 0(2 3 • 5 3 • 3 2 ) = (2 3 - 2 2 )(5 3 - 5 2 )(3 2 - 3 1 ) = 
4 • 100 • 6 = 2400 
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Chapter 19 

The Chinese Remainder 
Theorem 


19.1 Definition A linear congruence has the form ax = b (mod n) where x is 
a variable. 

19.2 Example The linear congruence 2x = 1 (mod 3) is solved by x = 
2 because 2-2 = 4 = 1 (mod 3). The solution set of that congruence is 

{...,2,5,8,11,-..}. 

19.3 Example The congruence 4x = 1 (mod 2) has no solution, because 4x 
is even, and so is not congruent to 1, modulo 2. 

19.4 Lemma Fix a modulus m and a number a. The congruence ax = b 
(mod m) has a solution if an only if gcd(a, in) \ b. If a solution x 0 does exist 
then, where d = gcd(a, b), the set of solutions is 


{.. •, x 0 + (—m/d),Xo, x 0 + ( m/d),x 0 + (2 m/d),x 0 + (3 m/d ),... } 


the residue class [xo] modulo m/d. 

PROOF. The existence of an x solving ax = b (mod m) is equivalent to the 
existence of a k such that ax — b = km, which in turn is equivalent to the 
equivalence of a fc such that xa + (— k)m = b. With that, this result is a 
restatement of Lemma 5.7. QED 

One generalization of Lemma 19.4 is to consider systems of linear congru¬ 
ences. In 1247, Ch’in Chiu-Shao published a solution for a special case of that 
problem. We first need a preliminary result. 

19.5 Lemma If gcd(a, b) = 1 and c is a number such that a \ c and b \ c then 
ab | c 
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PROOF. Because a \ c and b \ c there are numbers k a ,kb such that k a a = c and 
kbb = c. By Bezout’s Lemma, there are s and t such that as + bt, = 1. Multiply 
by c to get cas + cbt = c. Substitution gives ( kbb)as + (k a a)bt = c. Then ab 
divides the left side of the equation and so ab must divide the right side, c. QED 

19.6 Theorem (Chinese Remainder Theorem) Suppose that mi,. .. ,m n 

are pairwise relatively prime (that is, gcd(m.j,my) = 1 whenever i ^ j). Then 
the system of congruences 


x = a\ (mod mi) 
x = a 2 (mod m 2 ) 


x = a n (mod m n ) 

has a unique solution modulo ?i 7 .i?n 2 ... m n . 

proof. Let M = ... m n and for i G {l,...,n} let Mi = M/mi = 

miW 2 ... ... m n . Observe that gcd(Afi, ?n,;) = 1 and so Lemma 19.4 

says that the linear congruence ARx = 1 (mod mi) has a set of solutions that 
is a single congruence class [xi] modulo m,j. 

Now consider the number 


s 0 = aiMiaq + a 2 M 2 x 2 + • • • + a n M n x n . 

We claim that Sq solves the system. For, consider the i-th congruence x = a.i 
(mod ?n,;). Because mi divides Mj when i ^ j, we have that So = OjMjXj 
(mod mi). Since Xi was chosen because of the property that M.;X, = 1 (mod mi), 
we have that so = at ■ 1 = a,: (mod m.j), as claimed. 

To finish we must show that the solution is unique modulo M. Suppose that 
x also solves the system, so that for each i G { 1 ,..., n} we have that x = ai = x 0 
(mod mi). Restated, for each i we have that n* | ( x — Xq). 

We can now show that miro 2 .. .m n \ (x — xq). We have that gcd(mi,m 2 ) = 
1 and mi \ (x — a;o) and m 2 | (x — xq), so the prior lemma applies and we 
conclude that mim 2 | (x — Xo). In this way, we can build up to the entire 
product m-i... m n . QED 



Chapter 20 

Fermat’s Little Theorem 


20.1 Definition For [a] £ U m , the powers of the residue class are given by 
[a ] 1 = [a], [a ] 2 = [a][a], etc. 

20.2 Lemma If [a] £ U m then [a] n € U m for n > 1, and [a]" = [a"]. 

PROOF. We will check this by induction on n. The n = 1 base case is trivial: 
[a ] 1 = [a] = [a 1 ], and by assumption [a] € U m . For the inductive step, suppose 
that [a] k = [a fe ] G U m for k > 1 and consider the k + 1-st power. 

[«] fc+1 = H fe [a] = [a fc ] [a] = [a k a] = [a k+1 } 

By induction the theorem holds for all > 1. QED 

20.3 Theorem (Euler’s Theorem) If m > 0, and a is relatively prime to to, 

then = 1 (mod to). 

PROOF. For to > 0, we have that gcd(a, to) = 1 if and only if [a] £ U m . The prior 
result gives that a n = 1 (mod to) [ a n } = [1] [ a] n = [1]. Therefore, 

Euler’s Theorem is equivalent to the following: if to > 0 and [a] £ U m then 

[ a ]0(m) = [i]. 

We will write X- t . X- 2 , ■ ■ ■. X ^( m ) for the residue classes in U m . 

We first show that if X £ U m then the set O = {XXi, XX 2 , ■ ■ ■ ,XX^ m ^} 
equals the set U m . Containment one way is easy: any member of O is a member 
of U m by the closure property of Theorem 18.7. For containment the other 
way, consider Xi £ U m ,and note that Theorem 18.9 shows that the equation 
X 0 x = Xi has a solution x = Xj for some j, so X, = XXj is an element of O. 

Next, for any X £ U m consider the product XXiXX 2 ■ ■ ■ XX^^. The 
associative property says that we can parenthesize this term in any way, and 
the prior paragraph then gives that the product is {XXi){XX 2 ) ■ ■ ■ ( XX ^( m )) = 
X±X 2 • • • Xf^^jnp 

Finally, let A = XiX 2 ■ ■ ■ and for any X £ U m consider X^^A. 

The commutative property of Theorem 18.7 gives that 

X<P(m) A = X^X x X 2 ■ ■ ■ X H m) = (XXJiXXi) ■ ■ ■ (XX Hm) ). 
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The prior paragraph then shows that X^ m ^A = A. 

Multiplying both sides of that equation by the inverse A* of A gives 

^X^ m ^A)A* = X^ m \AA*) = X^ m )[l] = X 

on the left and AA* = [1] on the right, as desired. QED 

20.4 Example Fix m = 12. The positive integers a < m with gcd (a,m) = 1 
are 1, 5, 7 and 11, and so </>(m) = 4. We will check Euler’s result for all four. 

First, l 4 = 1 (mod 12) is clear. Next, 5 2 = 1 (mod 12) since 12 | 25 — 1, and 
so 5 4 = (5 2 ) 2 = l 2 (mod 12). From that one, and because 7 = —5 (mod 12) 
and 4 is even, 7 4 = 5 4 (mod 12) = 1 (mod 12). And, fourth, 11 = —1 (mod 12) 
and again since 4 is even we have that ll 4 = (—l) 4 (mod 12) = 1 (mod 12). 

20.5 Theorem (Fermat’s Little Theorem) If p is prime, and a is relatively 
prime to p, then a v ~ 1 = 1 (mod p). 

PROOF. Where p is prime, <(>(p ) = p — 1. QED 

20.6 Example Fermat’s Little Theorem can simplify the computation of 
a" mod p where p is prime. Recall that if a n = r (mod p) where 0 < r < p, then 
a n mod p = r. We can do two things to simplify the computation: (i) replace a 
by a mod p, and (ii) replace n by n mod (p — 1). 

Suppose that we want to calculate 1234 7865435 mod 11 Note that 1234 = 
—1 + 2 — 3 + 4 (mod 11), that is, 1234 = 2 (mod 11). Since gcd(2,11) = 1 we 
have that 2 10 = 1 (mod 11). Now 7865435 = (786543) • 10 + 5 so 

2 7865435 = 2 (786543).10+5 (mod n) 

= ( 2l0 ) 78 6543. 2 5 ( modn ) 

= l 786543 . 2 5 (mod 11) 

= 2 5 (mod 11), 

and 2 5 = 32 = 10 (mod 11). Hence, 1234 7865435 = 10 (mod 11). It follows that 
1234 7865435 mod 11 = 10. 

20.7 Remark Fermat’s theorem is called “little ” as a contrast with Fermat’s 
Last Theorem, which states that x n + y n = z n has no solutions x,y,z € N 
when n > 2. For many years this was the most famous unsolved problem in 
Mathematics, until it was proved by Andrew Wiles in 1995, over 350 years after 
it was first mentioned by Fermat. Fermat’s Little Theorem is much easier to 
prove, but has more far-reaching consequences for applications to cryptography 
and secure transmission of data on the Internet. 



Chapter 21 

Probabilistic Primality 
Tests 


Fermat’s Little Theorem says that if p is prime and 1 < a < p — 1, then a p 1 = 1 
(mod p). It has this converse. 

21.1 Theorem If rn > 2 and for all a such that 1 < 8 < m - 1 we have 
a m ~ i = 1 (mod to) then to must be prime. 

PROOF. If the hypothesis holds, then for all a with 1 < a < m — 1, we know 
that a has an inverse modulo to, namely, a m ~ 2 . By Theorem 18.3, this says 
that for all 1 < a < to — 1 we have that gcd(a, in ) = 1. But this means that in 
is prime, because if not then we would have to = ab with 1 < a, b < to, which 
would mean gcd(a, in) = a > 1. QED 


Therefore, one way to check that a number to is prime would be to check that 
l m— i = 1 (mod to), and that 2 m_1 = 1 (mod to), ..., and that in — l m_1 = 1 
(mod in ). 

This check is a lot of work, but it does have an advantage. Consider to = 63. 
Note that 2 6 = 64 = 1 (mod 63) and raising both sides to the 10-th power gives 
2 60 = 1 (mod 63). Multiplying both sides by 2 2 yields the conclusion that 
2 62 = 4 (mod 63). Since 4^1 (mod 63) we have that 2 62 ^ 1 (mod 63). This 
tells us, without factoring 63, that 63 is not prime. 

On the other hand, knowing only that 2 m ~ 1 = 1 (mod to) is not enough 
to show that in is prime. For instance, 2 m_1 = 1 (mod to) for the composite 
number to = 341. 

Nonetheless, consider only the base 6 = 2. There are 455,052,511 odd primes 
p < 10 10 , all of which satisfy 2 P ~ 1 = 1 (mod p). There are only 14,884 com¬ 
posite numbers 2 < m < 10 10 that satisfy 2 m ~ 1 = 1 (mod to). Thus, for a 
randome number m with 2 < to < 10 , if to satisfies 2 m_1 = 1 (mod in) then 
the probability that m is prime is 


455,052,511 
455,052,511 + 14,884 


.999967292. 
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In other words, if we find that 2 m_1 = 1 (mod in), then it is highly likely (but 
not a certainty) that m is prime, at least when m < 10 10 . Thus we are led to 
the following algorithm (expressed in the syntax of Maple). 

> is_prob_prime:=proc(n) 

if n <=1 or Power(2,n-l) mod n <> 1 then 
return "not prime"; 
else 

return "probably prime"; 
end if; 
end proc: 

What happens if we use 3 instead of 2 in the above probabilistic primality 
test? Or, better yet, what if we evaluate a m_1 mod in for several different a’s? 

The number of primes less than 10 6 is 78,498. The number of numbers 
to < 10 6 that are composite and such that 2 m_1 = 1 (mod to) is 245. The 
number of numbers to < 10 6 that are composite and such that both 2 m ~ 1 = 1 
(mod in) and = 1 (mod to) is 66 . The number of numbers in < 10 6 that 
are composite and such that a m ~ 1 = 1 (mod to) where a is any of the first 
thirteen primes is 0. (If to > 10 6 and a m_1 = 1 (mod to) for all a in the set of 
the first thirteen prime then it is highly likely, but not certain, that to is prime.) 

That is, if we check for primality by using the scheme of this chapter then we 
may possibly find out early that the number is not prime, having done very little 
work. Otherwise, as we work our way through bases a £ [1..to), calculating 
whether a m ^ 1 = 1 (mod to), we gain confidence that to is prime. This is the 
Solovay-Strassen pseudoprimality test. 

In practice, there are better probabilistic primality tests than the one de¬ 
scribed here. For instance, the built-in Maple procedure isprime is a very 
sophisticated probabilistic primality test. So far no one has found an integer n 
for which isprime(n) gives the wrong answer. 



Chapter 22 

Representations in Other 
Bases 


22.1 Definition Let b > 2 and n > 0. The base b representation of n is n = 
[a k ,a k - i,... ,ai,ao] b for some k > 0, where n = a k b k + ak-ib 1 *- 1 + ■ • - + ai&+ao 
and a, G {0, 1 ,..., b — 1 } for i = 0, 1 ,..., k. 

22.2 Example (1) 267 = [5,3,1] 7 , since 267 = 5 • 7 2 + 3 • 7 + 1 

(2) 147 = [1,0,0,1,0,0,1,1] 2 , since 147 = 1 • 2 7 + 0 • 2 6 + 0 • 2 5 + 1 • 2 4 + 0 • 
2 3 + 0 - 2 2 + 1 - 2 + 1 

(3) 4879 = [4, 8 ,7,9]i 0 , since 4879 = 4 • 10 3 + 8 • 10 2 + 7 • 10 + 9 

(4) 10705679 = [A, 3,5, B, 0, F]i 6 , since 10705679 = 10 • 16 5 + 3 • 16 4 + 5 • 
16 3 + 11 • 16 2 + 0-16 + 15 

Observe that a number’s base 10 representation is just its ordinary one. 
The representations are said to be in binary if b = 2, in ternary if b = 3, 
in octal if b = 8 , in decimal if b = 10, and in hexadecimal if b = 16. If b is 
understood, especially if b = 10 , we write a k a k -± • • • a iOo, without the subscript 
base. In the case of b = 16, which is used frequently in computer science, for 
the ai of 10, 11, 12, 13, 14 and 15 we use A, B , C, D, E and F, respectively. 

For a fixed base b > 2, the numbers s the digits of the base b representa¬ 
tion. In the binary case, the a,;’s are bits, a shortening of “binary digits”. 

22.3 Theorem If b > 2 then every n > 0 has a unique base b representation. 
PROOF. To show that a representation exists, iterate the Division Algorithm: 

n = bq 0 + r 0 0 < r 0 < b 

qo = bqi + r 1 0 < ri < b 

qi = bq 2 + r 2 0 < r 2 < b 

q k = bq k +\ + rfc+i 0 < r k+ 1 < b. 
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Note that n > qo > qi > ■ ■ ■ > qu- This shows that iteration of the Division 
Algorithm cannot go on forever, and we must eventually obtain qe = 0 for 
some £, so that qe-i = b • 0 + r+ We claim that the desired representation is 
n = [re,re-i, ... ,ro]. For, note that n = bqo + ro and qo = bq± + r\, and hence 
n = b{bqi + r i) + ?’ 0 = b 2 qi + br\ + tq- Continuing in this way we find that 
n = b t+1 qe + b^rg + • • • + br\ + ro- And, since = 0 we have 

(*) n = b e r£ + ■ ■ ■ + br 1 + r 0 , 

which shows that n = [rg,... ,ri,ro] b - 

To see that this representation is unique, note that from equation (*) we 
have 

n = b (b l ~ 1 rt H- 1 - ri) + r 0 , 0 < r 0 < b. 

Because ro is uniquely determined by n, so is the quotient q = b t ~ 1 rc H-b ri. 

A similar argument shows that T\ is uniquely determined. Continuing in this 
way we see that all the digits re, re- 1 ,..., ro are uniquely determined. QED 

22.4 Example We find the base 7 representation of 1,749. 

1749 = 249-7 + 6 
249 = 35 • 7 + 4 
35 = 5 • 7 + 0 
5 = 0 • 7 + 5 


Hence 1749 = [5,0,4, 6 ] r . 

22.5 Example This finds the binary representation of 137. 

137 = 2 • 68 + 1 
68 = 2 • 34 + 0 
34 = 2 • 17 + 0 
17 = 2-8 + 1 
8 = 2 • 4 + 0 
4 = 2 • 2 + 0 
2 = 2 • 1 + 0 
1 = 2 • 0 + 1 


Therefore 137 = [1,0, 0, 0,1,0, 0,1] 2 . 

22.6 Remark We can sometimes “eyeball” the representation in another base 
of a small number. For instance, we can see how to represent n = 137 in binary, 
without the machinery of the proof. Note that 2 1 = 2, 2 2 = 4, 2 3 = 8 , 2 4 = 16, 
2 5 = 32, 2 6 = 64, 2 7 = 128, and 2 8 = 256. By eye, we spot that the value closest 
to 137 but not greater than it is 2 7 , and we compute that 137 — 2 7 = 9. The 
power of 2 closest to it but not above 9 is 2 3 , and 9 —2 3 = 1. Finally, 1 is a power 
of 2, since 1 = 2°. Therefore 137 = 2 7 + 2 3 + 2° and so 137 = [1,0, 0, 0,1,0,0,1] 2 . 



Chapter 23 


Computation of a 


N 


mod m 


Some Number Theory work involves computing with large numbers. Since com¬ 
puter multiplication of numbers is a slow operation (relative to computer addi¬ 
tion), we can ask: where n is any positive integer, what is the smallest number 
of multiplications required to compute a n l 

For instance, the naive way to calculate 2 8 is to do seven multiplications. 


2 2 = 2 • 2 = 4 

2 3 = 2 • 4 = 8 

2 4 = 2 • 8 = 16 

2 5 = 2 • 16 = 32 

2 6 = 2 • 32 = 64 

2 7 = 2 • 64 = 128 

2 8 = 2 • 128 = 256 

In general, computing a n by this naive method requires n — 1 multiplications. 
But we can compute 2 s with only three multiplications 


2 2 = 2 • 2 = 4 

2 4 = (2 2 ) 2 = 4-4 = 16 

2 8 = (2 4 ) 2 = 16 • 16 = 256 

If the exponent has the form = 2 fc then this successive squaring method 
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requires only fc-many multiplications. 


a = a ■ a 
, 2 2 _ ( n 1\1 


2 ( 2\2 2 2 
a = (a ) = a • a 


a 


= (a 2 r = 


This is quite a savings because if n = 2 k then k is generally much smaller than 
n — 1, just as 3 is smaller than 7. 

This is the foundation of the binary method to compute a n . It is best 
explained by example. 

23.1 Example To compute 3 15 , first express the exponent in binary 15 = 
2 3 + 2 2 + 2 + 1 = [1,1,1,1] 2 . Thus, 3 15 = 3 23 • 3 2 " • 3 2 • 3. 

Next, we use successive squaring to get the factors in that expansion of 3 15 . 

3 2 = 3 • 3 = 9 

3 22 = 9 • 9 = 81 

3 23 = 81 • 81 = 6561 
Putting those factors together 

3 • 3 2 = 3 • 9 = 27 

(3 • 3 2 ) • 3 22 = 27-81 = 2187 

(3 • 3 2 • 3 22 ) 3 23 = 21 87 • 6561 = 14348907 

gives that 3 15 = 14348907. This took just six multiplications, while the naive 
method would have taken fourteen. (Finding the binary representation of 15 
took some extra effort, but not much.) 

23.2 Theorem Computing x n using the binary method requires |_lg(n)J divi¬ 
sions and at most 2[lg(n)J multiplications. 

PROOF. If n = [a r ,..., ao ]2 and a r = 1 then 2 r < n << 2 r+1 . By the familiar 
properties of any logarithm, lg(2 r ) < lg(n) < lg(2 r+1 ). Since lg 2 (2 x ) = x this 
gives r < lg(n) < r + 1, hence r = |_lg(n.)J. Note that r is the number of times 
we need to divide to get n’s binary representation n = [a r ,..., 00 ) 2 - 

To compute the powers x, x 2 , x 2 ,..., x 2 by successive squaring requires 
r = |_lg( n )J multiplications and similarly to compute the product 

j,2 r . ^,a r -l2 r-1 , . , , j,ao 


requires r multiplications. So after obtaining the binary representation we need 
at most 2r = 2[lg(n)J multiplications. QED 
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Note that if we count an application of the Division Algorithm and a mul¬ 
tiplication as having the same cost then the above tells us that we need at 
most 3|_lg(n)J operations to compute x n . So, for example, if n = 10 6 , then 
3Llg(n)J=57. 

To compute a n mod m, we use the binary method of exponentiation, with 
the added refinement that after every multiplication we reduce modulo m. This 
keeps the products from getting too big for our computer or calculator. 

23.3 Example We compute 3 15 mod 10: 

3 2 = 3- 3 = 9 = 9 (mod 10) 

3 4 = 9 • 9 = 81 = 1 (mod 10) 

3 8 = 1-1 = 1 = 1 (mod 10) 

and so 3 15 = 3 8 • 3 4 • 3 2 • 3 1 = 1 • 1 • 9 • 3 = 27 = 7 (mod 10). 

Note that 3 15 = 7 (mod 10). In Example 23.1 we calculated that 3 15 = 
14348907 which is clearly congruent to 7 mod 10, but the multiplications there 
were not so easy. 

23.4 Example To find 2 644 mod 645, we first get the binary representation 
644 = [1,0,1,0, 0,0,0,1, 0,0] 2 . That is, 644 = 2 9 + 2 7 + 2 2 = 512 + 128 + 4. By 
successive squaring and reducing modulo 645 we get 

2 2 = 2- 2 = 4 = 4 (mod 645) 

2 4 = 4 • 4 = 16 = 16 (mod 645) 

2 8 = 16 • 16 = 256 = 256 (mod 645) 

2 16 = 256 • 256 = 65,536 = 391 (mod 645) 

2 32 = 391 • 391 = 152, 881 = 16 (mod 645) 

2 64 = 16 • 16 = 256 = 256 (mod 645) 

2 128 = 256 • 256 = 65,536 = 391 (mod 645) 

2 256 = 391 • 391 = 152, 881 = 16 (mod 645) 

2 512 = 16 • 16 = 256 = 256 (mod 645). 

Now 2 644 = 2 512 • 2 128 • 2 4 , and hence 2 644 = 256 • 391 • 16 (mod 645). So 
256 • 391 = 100099 = 121 (mod 645) and 121 • 16 = 1936 = 1 (mod 645). Hence 
2 644 mod 645 = 1. 
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Chapter 24 


Public Key Cryptosystems 


Everyone has tried secret codes. A common one is the Caesar cipher: the 
sender and the recipient agree in advance to express letters as numbers (1 for 
A, 2 for B , etc.) and also agree to use an encoding that offsets the message; 
for instance f(n) = (n + 13) mod 26 offsets the letters by 13. The sender then, 
in place of transmitting the number n, will transmit f(n )—instead of A, the 
sender will transmit K, the thirteenth letter. This code is very easy to break, 
but nonetheless notice that there is a general encryption/decryption scheme of 
sending offset letters, and within that scheme it relies on the single secret key, 
the 13. 

In 1976, W Diffie and M Heilman proposed a new kind of cryptographic 
system where there are two keys. A message encrypted with the first key can be 
decrypted with the second, and a message encrypted with the second key can 
be decrypted with the first. We will first illustrate some advantages of such a 
system and then give one way to produce such key pairs. 

24.1 Example If two people, Alice and Bob, want to have private commu¬ 
nications then each can take a key. Bob alone can read Alice’s messages, and 
Alice alone can read Bob’s. 

24.2 Example Alice can keep one key a secret, and publish the other key in a 
public place such as the Internet. Then people who receive an encrypted message 
that claims to be from Alice can get Alice’s public key and try to decrypt the 
message. If the result is sensible text, then Alice must have been the one who 
encrypted it, since she kept her other key private. This is authentication ; her 
message has been digitally signed. 

Also, people who want to send a message to Alice in private can encrypt it 
with her public key. Only she can decrypt it, using her private key. 

24.3 Example Key pairs can be used to do things that seem impossible. 
Suppose that Alice and Bob want to settle a dispute by flipping a coin, but they 
must do so over the Internet. Each person will flip separately, and they agree 
that Alice wins if the two coins come out the same while Bob wins if they are 
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different. However they do not trust each other and so they cannot just email 
each other the results. How can they agree if neither will believe the other? 

Each person generates a key pair. Each then sends the other the message 
with “heads” or “tails” encrypted using one of their two keys. After that, each 
person publishes their other key, the one that they did not use to encrypt. The 
other person can now decrypt the message they’ve received — they are sure that 
they are not being cheated because they now have the other person’s outcome, 
albiet encrypted (and the key pairs have the property that finding a new key 
pair that makes the message decrypt the other way is essentially impossible). 

Implicit in these examples are a number of technical requirements on key 
pairs: from either key we should not be able to find the other, we should not be 
able to decrypt the message by just trying every possible key, etc. These techni¬ 
cal requirements have been met by a number of schemes. The most important 
is RSA, due to R Rivest, A Shamir, and L Adelman in 1977 [11]. This chapter 
outlines its number-theoretic underpinning. 

Assume that our message has been converted to an integer in the set J m = 
{0,1,2,... , m — 1} where m is some positive integer to be determined. (For 
example, we can take the file as a collection of bits and interpret it as a num¬ 
ber written in binary.) Generally this is a large integer. We will require two 
functions: 

E : J m —> J m ( E for encipher) and D : J m —> J m ( D for decipher). 

By ‘encipher’ and ‘decipher’ we only mean that D(E{x)) = x for all x £ J m . 
We first need two statements about congruences. 

24.4 Lemma Let mi,m 2 £ Z+ be relatively prime. Then a = b (mod mi) 
and a = b (mod m 2 ) if and only if a = b (mod TO1TO2). 

PROOF. One direction is easy: if a = b (mod mi m2) then there is k £ Z such 
that a — b = k{mim< 2 ). Rewriting that as a — b = (fcmi)m2 shows that a — b is 
a multiple of m2 and so a = b (mod m2). The other equivalence is similar. 

If a = b (mod mi) and a = b (mod m 2 ) then there are ki,k 2 £ Z such 
that a — b = Aqmi and a — b = fem^- Therefore k\m\ = ^2^2- This shows 
that mi | As mi is relatively prime to m2, Lemma 5.6 then gives that 

mi | /c2- Writing fc 2 = km\ for some k, and substituting into the earlier equation 
a — b = /C2W2 gives that a — b = fcm.im.2. Therefore a — b \ miiri 2 and so a = b 
(mod mi m2). QED 

24.5 Lemma Let p and q be two distinct primes and let m = pq. Suppose 
that e and d are positive integers that are inverses of each other modulo < f>(m ). 
Then x ed = x (mod m) for all x. 

PROOF. By Theorem 18.17, <j>(m ) = (p — l)(g — 1). Since ed = 1 (mod < j>(m )) 
we have that ed — 1 = k(j>(m ) = k(p — l)(g — 1) for some k. Note the k > 0 
unless ed = 1, in which case the theorem is obvious. So we have 

(*) ed = kcj)(m) + 1 = k{p — l)(q — 1) + 1 
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for some k > 0. 

We will show that x ed = x (mod b) for all x. There are two cases. For 
the first case, if gcd(a;,p) = 1 then by Fermat’s Little Theorem we have that 
x p_1 = 1 (mod p). Raising both sides of the congruence to the power (q — 1 )k 
gives = 1 (mod p). Then multiplying by x gives x^ p ~ 1 ' >l ' q ~ 1 ' ,k+1 = x 

(mod p). That is, by (*) 

(**) x ed = x (mod p). 

For the second case, the gcd(cc,p) = p case, the relation (**) is obvious, since 
then x = 0 (mod p). 

A similar argument proves that x ed = x (mod q) for all x. So by Lemma 24.4 
and the fact that gcd(p, q) = 1, we have that x ed = x (mod to) for all x. QED 

24.6 Theorem Let p and q be two distinct primes, let to. = pq , and suppose 
that e and d are positive integers that are inverses of each other modulo <j){rn). 
Where J m = {0,1,2,..., to. — 1}, define E : J m —> J m and D : J m —> J m by 

E{x) = x e mod to. and D(x) = x d mod to. 

Then E and D are inverse functions. 

PROOF. It suffices to show that D(E(x)) = x for all x £ J m . Suppose that 
x £ J m , and that E(x) = x e mod m = j*i, and also that D (ri) = r d mod to = 
T 2 - We must show that r2 = x. Since x e mod to = r\ we know that x e = r\ 
(mod m). Hence x ed = rf (mod to). We also know that r d = r 2 (mod m) and 
hence x ed = r 2 (mod to). By Lemma 24.5, x ed = x (mod to) so we have that 
x = r 2 (mod to). Since both x and r 2 are elements of J m , both are the principle 
residue: x = r 2 - QED 
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Appendix A 

Proof by Induction 


Most of the proof methods used in mathematics are instinctive to a person with 
a talent for the work. This section covers a method, the method of Mathematical 
Induction that is not. 

As with all proofs, we will have some assertion to prove. Each assertion will 
say that something is true for all integers. Thus, we can denote the assertion 
P(ri). Our first example is the proof that for all n, if n > 5 then 2 n > 5n. 

P(n): n > 5 <t=> 2" > 5n 

An argument by induction involves two steps. In the base step we show that 
P is true for some first integer. Typically, that is a straightforward verification. 
For our example, we show that P(5) is true by just checking that 2 5 = 32 is 
indeed greater than 5 • 5 = 25, which of course it is. 

The second step is called the inductive step. We must show that 

if P{ 5), ..., P(k ) are all true then P(k + 1) is also true. 

At the end of the proof we will show why this suffices. For the moment note 
only that we are not asserting that P(5), ..., P{k) are in fact all true (as that 
would be assuming the thing that we are to prove); instead we are proving that 
if they are true then P(k + 1) follows. 

To prove this if-then statement, take the inductive hypothesis that P(5), 
..., P(k) hold. Then, by the hypothesis that P(fc) is true we have 2 k > 5k, 
and Multiplying both sides by 2 gives 2 k+1 > lOfc We are trying to prove that 
2 k+1 > 5 (k + 1) so if we can show lOfc > 5k + 5 then we will be done. Because 
k > 5, we have that 5k > 5 and therefore lOfc = 5k + 5k > 5k + 5 = 5 (k + 1). 
We have therefore established P(k + 1) follows from the inductive hypothesis, 
as 2 k+1 > lOfc > 5{k + 1). That ends the inductive step. 

To see why the two steps togehter prove the assertion, note that we have 
checked the statement for 5. To see it is true for 6, note that in the inductive 
step we proved that P(5) <t=> P(6). To see that the statement is true for 7, note 
that we have proved in the inductive step that P(5) and P(6) P(7) (and the 
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prior sentence shows that P(6) holds). In this way we can see that the statement 
is true for all numbers n > 5. 

Here is an induction proof that is more streamlined, more like the ones given 
elsewhere in the book.. 

1.1 Proposition If n > 5 then 2" > 5 n. 

PROOF. We prove the proposition by induction on the variable n. 

If n = 5 then we have 2 5 > 5 • 5 or 32 > 25, which is true. 

Next, assume the hypothesis that 2 n > 5 n for 5 < n < k. Taking n = k gives 
that 2 fc > 5 k. Multiplying both sides by 2 gives 2 fc+1 > 10/c. Now 10/c = 5/c + 5/c 
and k > 5, and so 5 k > 5. Hence 10/c = 5/c + 5/c > 5/c + 5 = 5 (k + 1). It follows 
that 2 k+1 > 10/c > 5 (k + 1) and therefore 2 fc+1 > 5 (k + 1). 

Hence by mathematical induction we conclude that 2" > 5 n for n > 5. QED 



Appendix B 

Axioms for Z 


The set of natural numbers is N = {0,1, 2, 3, • • • }. The set of integers includes 
the natural numbers and the negative integers Z = 2,—1,0,1,2, •••}• We 

sometimes want to restrict our attention to the positive integers Z + = {1,2, • • • }• 

The rational numbers include all of the fractions Q = {n/m \ n,m £ 
Z and m/0}. The real numbers K. enlarge that set with the irrationals (which 
are too hard to precisely describe here). Note that Z + cNcZcQcR. 

In the first chapter we rely on some particularly important properties of Z, 
the axioms. 

1. If a, b £ Z, then a + b, a — b and ab £ Z. (That is, Z is closed under 
addition, subtraction and multiplication.) 

2. If a £ Z then there is no x £ Z such that a < x < a + 1. 

3. If a, b £ Z and ab = 1, then either a = 6= lora = 6= —1. 

4. Laws of Exponents For n, m £ N and a,b £ M. with a and b not both 0 we 
have (a n ) m = a nm and ( ab) n = a n b n . and a n a m = a n+m . 

5. Properties of Inequalities: For a, b, c in R. the following hold: if a < b and 
b < c, then a < c, and if a < b then a + c < b + c, and if a < b and 0 < c 
then ac < be, and if a < b and c < 0 then be < ac , and finally, given a and 
b, one and only one oi a = b , a < b, b < a holds. 

6. Well-Ordering Property Every non-empty subset of N contains a least 
element. 

7. Mathematical Induction Let P(n) be a statement concerning the integer 
variable n. Let no be any fixed integer. Then P(n) is true for all integers 
n > no if both of the following statements hold: (the base step) P(n) 
is true for n = no, and (the inductive step) whenever P(n) is true for 
n 0 < n < k then P(n) is true for n = k + 1. 
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Appendix C 

Some Properties of R 


3.1 Definition Where i£l, the floor (or greatest integer) [xj is the largest 
integer less than or equal to x. Its ceiling \x~\ is the least integer greater than 
or equal to x. 

For example, [3.1] = 3 and [3.1] = 4, [3J = 3 and [3] = 3, and [—3.1J = —4 
and [—3.1] = —3. 

From that definition we immediately have that [xj = max{?r € Z | n < x}, 
and that n = [xj •<=> n < x < n + 1. From this we have also that [x] < x and 
that [x] = x x€Z. 

3.2 Lemma (Floor Lemma) Where x is real, x — 1 < [xj < x. 

PROOF. Let n = [xj. Then by the above comments, we have n < x < n + 1. 
This gives immediately that [xj < x, as already noted above. It also gives that 
x < n + 1 which implies that x — 1 < n, that is, that x — 1 < [xj. QED 

3.3 Definition The decimal representation of a positive integer a is given by 
a = a n -ia n _2 ■ • • aicto where 

a = a„_ilO " _1 + a„_ 2 10™ -2 + • • • + ailO + ao 

and the digits a„_i, a n _ 2 , • ■ •, ai, ao are in the set {0,1, 2, 3,4,5, 6 , 7, 8 ,9}, with 
a n _i 7 ^ 0. This representation shows that a is, with respect to base 10, an n 
digit number (or is n digits long). 
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